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ABSTRACT 


Grain-size sampling of densely packed particulate materials 
by customary bulk sieve analysis is sometimes very difficult. In an 
attempt to find an equivalent alternative to bulk sieve analysis this 
thesis examines a probabalistic grain-size sampling theory proposed by 
Kellerhals et al (1975) here designated as the numerical method, and 
four related empirical sampling experiments. 

Three of the experiments incorporate results found within the 
literature. Data from Kellerhals' et al (1975) and Friedman's (1962) 
grid-by-number thin section experiments and McGinn's (1971) grid-by- 
number gravel bar investigation all tend to confirm the applicability of 
the numerical method and grid-by-number as a solution to this sampling 
problem. 

The fourth experiment involved the grid-by-number investigation 
of terrace gravels. It is shown that a terrace gravel surface may in 
some circumstances be treated as a thin section surface. The statistical 
results of this final study strongly confirm numerical method predictions. 

Generally, the results of these experiments indicate that for 
a wide variety of sampling situations the grid-by-number sampling tech- 
nique in conjunction with the numerical method can provide at least 


median and mean values equivalent to those obtained by sieving. 
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INTRODUCTION 


Grain-size sampling of densely packed particulate materials 
such as clastic sediments has usually depended upon customary bulk sieve 
analysis. In a number of sampling situations the application of this 
Standard technique is difficult. For instance, the disaggregation of a 
sandstone for sieve analysis or the collection of the large gravel 
samples that are necessary for representativeness may be so laborious as 
to be impractical. What is required and what has long been sought is an 
alternative system of sampling which avoids these impediments and yields 
grain-size distributions which are the same as those derived by bulk 
sieve analysis. 

Central to this thesis is the empirical investigation of one 
grain-size sampling theory which possibly provides an equivalent 
alternative to bulk sieve analysis. The probabalistic sampling theory 
proposed by Kellerhals et al (1975), here designated as the numerical 
method, employs a computer program to simulate grid-by-number sampling 
of a randomly selected thin section obtained from an isotropic material 
composed of identical ellipsoidal grains. For a given uniform material 
the numerical method predicts the relationships among true axial 
dimensions, square mesh sieve size and k values of the constituent 
equi-sized ellipsoidal grains. Furthermore, the probabalistic association 
of these ellipsoidal properties and the mean values of the major and 
minor apparent axial cumulative distributions are given. 

Mathematical models of complex natural structures and processes 
are generally idealizations. The numerical method assumes that the 


material being sampled is uniform, while virtually all clastic sediments 
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are nonuniform. For the numerical method to be useful it must be shown 
that its predictions are confirmed for the sampling of nonuniform materials. 
A number of sampling experiments employing nonuniform materials can be 
designed to test specific numerical method relationships. Fortunately, 
experimental data found within the literature can be utilized in some 
instances for this purpose. The thin section experiment accompanying 
the Kellerhals et al (1975) paper is reviewed and further analyzed. 
Data from Friedman's (1962) thin section and McGinn's (1971) gravel bar 
experiments are reexamined. In the three experiments numerical method 
predictions are tested through the comparison of grid-by-number and 
Sieve distributions. 

Other possibilities remain for checking numerical method pre- 
dictions. This study presents the results of a fourth experiment which 
employs terrace gravels. These gravels are often exposed in out-banks 
in which the face of the deposit is almost vertical. Proving that a 
thin section surface and a terrace gravel photograph are equivalent is 
essential for this sampling experiment and as such is the subject of an 
extensive discussion. The terrace gravel experiment involves the grid- 
by-number sampling and subsequent measurement of true axes in the field 
and apparent axes from SROLOaanNe. Or the same clasts. The numerical 
method predictions based on the true and apparent axial cumulative 


distributions are then examined. 
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LIST OF TERMS 


MPD denotes the minimizing plane dimension which is the line 
resulting from the intersection of the minimizing plane 
and the plane produced by sectioning. 

ASD denotes the apparent surface dimension which corresponds 
to the minor apparent axis of the grain's outline trace 
On the apparent surface 

Cag: denotes the center of gravity 

Y denotes the angle in degrees between the major axis and 
the line of sight passing through the center of gravity 

OGS denotes the overlapping grain surface which consists of 
an assemblage of surface grains belonging to all three 
contact modes. A familiar example is the surface of 
exposed gravel bars. 

EGS denotes the embedded grain surface which consists of an 
assemblage of surface grains belonging to contact modes 
2 and 3. A familiar example is the near vertical face 
associated with terrace gravels. 

A denotes the center of gravity - surface plane distance 


up to and including one radius 


n denotes a number or size of unspecified value or magnitude 
D denotes the linear size of cubic grains 

V _ denotes the specimen volume 

S denotes a smooth planar specimen surface 

A denotes the area of S 
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sieve 


Group 1 


Group 1, 


denotes the total area of S covered by cubic grains with 
linear size Dd. 

denotes the combined lengths of cubic grains with linear 
size D, which are touched by a transect placed on S 
denotes the number of cubic grains with linear size D,, 
which fall under the intersections of a grid placed on S 
denotes the maximum dimension of a square hole of side ~ 
length D.. This dimension is along the two diagonals 
whose length Ly = 206 
Subscript t signifies that the term's value depends upon 
the measurement of surface grain true axes. 

Subscript p signifies that the term's value depends upon 
both experimental data and numerical method predictions. 
Subscript sieve signifies that the term's value depends 
upon the results of customary bulk sieve analysis. In 
the case of quadrant-by-weight the sieve subscript denotes * 


that this distribution has been converted to its volume- 


by-weight equivalent. 


True Axial Terms 
A,B, and C denote the major (large), intermediate, and 
minor (small) axes respectively of an ellipsoid. 
A.s Bi. 
A, B, and C axes. 


C, denote the field measurement of a grain's true 


AB, > ABC, » and BC, denote different arithmetic means 
derived from a grain's A.» B. and C (ABC, is the triaxial 
mean of a grain). 
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A AB B ABC BC 0° and Cis denote various 


t50° “2yt50.. -t50" LOO aro 0 
medians of grid-by-number distributions based on field 


measurements of true axes. 


ieee, 


A. AB,» B,» ABC,, BC,, and C, denote various means of 


grid-by-number distributions based on field measurements 


of true axes. 


Group 1 C50 and Es denote the predicted median and mean minor 


axis length whose values are based on apparent axial 


values and numerical method predictions. 


b 
Example computation of C_,_,: Using k, = £90. in Figure 14b), 
p50 2 a+50 


(B-C)100 _ e 
a = Vs Cogg = (1.0-¥) (disp) 


same way except Db, and a, are used. 


the mean counterpart of C can be calculated in the 


Apparent Axial Terms 


Group 2 a and b denote the major and minor axes of an ellipse or 
the major and minor apparent axes of an ellipsoid. 

a and b denote the means of apparent axial distributions 
based on a and b. 

Group 2+ ay and b, denote the major and minor apparent axes 
respectively of a surface grain associated with either 
an apparent or thin section surface. 
ab, denotes the arithmetic mean derived from a, and by. 
A450? ab. cg» and bi 5g denote various medians of grid-by- 


number distributions based on measurements of surface 


grain apparent axes. 
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Group 


Group 


Group 3 


2 


3 


sieve 


ay: ab,. and b, denote various means of grid-by-number 
distributions based on measurements of surface grain 
apparent axes. 

4050 and ay denote the predicted median and mean apparent 


major axis length whose values are based on true axial 


values and numerical method predictions. 


B50 C 


: j t50 
Example computation of a_,,: Using k,=,—— and k,=;—— 
p50 ire 2 Bicy 
ap a-B)100 _ y. . wh 
in Figure 14a), arse ven 450 Beco (T69 1) 


ay the mean counterpart of 250 can be calculated in the 
same way except Ay» By and C, are used. 

bo50 and bo denote the predicted median and mean minor 
apparent axis length whose values are based on true axial 


values and numerical method predictions. 


C 
Example computation of b50° Using ko=gee in Figure 14b), 
t50 


(b-C)100 _ Y Wy, 


Cate’ Onso1a 50 100,” 


p50 


bo the mean counterpart of b550 can be calculated in the 


same way except B, and C, are used. 


Sieve Terms 
an denotes the square mesh sieve diameter (the side length 
of a square hole). 


D awd denotes the volume-by-weight sampling procedure 


(customary bulk sieve analysis). 


D denotes the median of a grain-size distribution 


sieve50 
produced by bulk sieve analysis. 


Dee denotes the mean of a grain-size distribution 


produced by bulk sieve analysis 
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Group 3 


o 


D denotes the bulk sieve analysis of the coarse 


sieve2 
(> 8mm.) portion of the subsurface sample. 


De ave D denotes the mean of the grain-size distribution 


obtained from Deeye?: 
Q denotes the quadrant (area)-by-weight sampling procedure 
whereby all surface grains (> 8mm.) within a specified 
area are removed and sieved. 

Q denotes the mean of the grain-size distribution obtained 
from Q. 


Q 
sion of Q to volume-by-weight by the weighting factor a 


Seve denotes the distribution resulting from the conver- 


denotes the mean of the Q grain-size distribution. 


sieve sieve 


4550 and de denote the predicted median and mean sieve 
diameter whose values are based on apparent axial values 


and numerical method predictions. 


¢ 
; ; _ “p50 2\5 
Example computation of q950° dsp572k, 2(1+k2)] ’ 


where kz 

— t50 ee oe ah 
d D is calculated in the same way except by» ays and C, 
are used. 

D 


sp50 and DS denote the predicted median and mean sieve 


diameter whose values are based on true axial values and 
numerical method predictions. 


: : t50 2\5 % 

Example corte ton of Deoso! 5950°2k> [2(1+k5)] * 
ee OU 
where ko-3—— 
t50 


Dep is calculated in the same way except B, and C, are used. 
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CHAPTER ONE 


SURFACE SAMPLING EQUIVALENCE PROBLEMS 


1.1 General Introduction 


Central to this thesis is the examination of the Kellerhals 


et al(1975) numerical method and four grain-size sampling experiments 


which test its predictions. This introductory chapter discusses four 


topics of critical importance to the numerical method and its empirical 


tests: 


1) 


2) 


3) 


The numerical method assumes a densely packed granular 
material. It is essential that properties of this material 
relevant to grain-size sampling be identified and defined. 
The numerical method makes grain-size predictions specifically 
for thin section surfaces. While two of the experiments 
utilize thin sections the other two use gravel surfaces, one 
associated with the surface layer of exposed gravel bars, 

the other with terrace gravels. It is necessary to evaluate 
the characteristics of these surfaces relative to that of a 
thin section. 

Sampling procedures which require only the surface of a 
deposit are quite different from the standard volumetric 
procedure, bulk sieve analysis. A surface sampling procedure 


should be geometrically equivalent to this standard. 
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4) In customary bulk sieve analysis the square mesh 
sieve size De» is the grain-size measure. The 
relationships between D. and true or apparent axial 
grain-size measures frequently used in surface 


sampling procedures are investigated. 


1.2 Properties of a Densely Packed Granular Material 
1.2.0 Introduction 


Granular materials may be described as being densely packed or 
dilutely distributed (Kellerhals et al 1975). Earth sciences are usually 
concerned with densely packed grains in contact, while dilutely distri- 
buted granular materials are of greater interest to biologists. Since 
only densely packed material is examined here, henceforth it is simply 
termed material. The following discussion defines those textural 
characteristics of a material relevant to grain-size analysis. 

12.1 “SOntING 

The constituent grains of an homogeneous material (Kellerhals 
and Bray, 1971, p.1175) may vary in size. Sorting is a measure of the 
degree of grain-size similarity. 

13202 prabuic 

Isotropic means having the same properties in all directions. 
As used here, a material possesses an isotropic fabric when the axes of 
constituent grains are ran¢ y oriented. Conversely, when there is a 


definite axial orientation the fabric is anisotropic. 
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1.2.3 Induration 

Induration is the process whereby a material is hardened into 
rock by exposure to heat, pressure or a cementing agent. Once indurated 
it becomes very difficult to disaggregate into constituent grains. 

1.2.4 Grain Shape 

Grain shape can be defined in terms of the ratio of the three 
axial dimensions. Using values for the large axis A, intermediate axis B, 
and the small axis C, the Zingg grain shape ratios ky apokom es can be 
obtained (Kellerhals et al 1975). Depending upon the values of k, and 
Ko a grain shape falls within one of four general classes a) tabular, 
oblate or discoidal b) equant, equiaxal or spherical c) bladed or 
triaxial d) prolate or rod-shaped (Whitten et al 1972). 

A wide range of sedimentary grains may be approximated by 
ellipsoids (Allen 1969; Kellerhals et al 1975). According to the latter 
reference the common range of k values for these ellipsoids is 0.55-0.75. 
Allen (1969) states that the shapes of these triaxial ellipsoids may be 
ppneoanened by k,=0.667, k5= 0.50. It is acknowledged that there is some 


discrepancy between these two sources, although this may be unimportant. 


1.3 Comparison of Gravel Bar, Terrace Gravel and Thin Section Surfaces 


L.3:0) Introduction 

Many textural properties of a material can be investigated by 
observing its surface. These properties include degree of sorting and 
homogeneity, presence of matrix, fabric and grain size. It is this latter 
property which is of concern here. 

In the following discussion, it is assumed that the material, 


the constituent grains, and the surface have certain basic properties. 
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Material 


The material is homogeneous, nonindurated and isotropic. 


Grains 


The constituent grains of this material are ellipsoidal in 
shape and have common k values 0.55-0.75. The grains may 


be of any size and the density of each grain mass is uniform. 


Surface 


The surface being sampled is approximately planar. The ideal 
plane associated with this surface is defined as the surface 
plane. If the surface is viewed from a given point on a 
line of sight normal to the surface plane, it is defined 

as the apparent surface. Grains observed on the apparent 
surface are termed surface grains. The apparent grain 
surface is defined as that part of the actual grain surface 
of a surface grain which can be viewed on the apparent 
surface. An apparent grain surface is bounded by its out- 
line trace. This is formed when an apparent grain surface 
appears to contact apparent grain surfaces of other adjacent 


surface grains. 


Employing these assumptions it is proposed that: 


1) 


there are three types of surfaces; overlapping grain surface 
(OGS), embedded grain surface (EGS) and thin section surface. 
the surface associated with the surface layer of gravel bars 
possesses an OGS and that of a near vertical face of a 
terrace gravel deposit possesses an EGS. 

grain-size analysis procedures developed for thin sections 


can be used on a terrace gravel EGS. 
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1.3.1 Inter Grain Contact 

The actual grain surface of a surface grain has one or more 
point contacts with at least one other grain. It is assumed that each 
contact point participates in supporting the grain in its present surface 
position. In well-sorted material the number of contact points tends to 
be low; conversely, if there is a matrix and the surface grain is 
relatively large, the number of contact points is probably much greater. 
A surface grain may be in contact with both matrix grains and similar 
sized grains. 

There are three basic classes of contact in which the surface 
grain may be involved; single point contact, two point contact, and three 
Or more point contact. In the second class a straight line may be 
envisaged to link the two points. In the third class any three contact 
points can flmetion as vertices of a triangular shaped figure. This 
figure forms a three point contact plane. 

The number of three point contact planes which can be formed 
from a set of n contact points can be determined by using C(n,r) which | 


is the number of r-subsets of a set of n elements. 


C(n,r) = nea fis where nsreNos me <angatly) 


) 
In this case r=3 Gitis3) = (123) 13! 


A three point contact plane may be extended until it is bounded 
by the actual grain surface. This plane is termed the extended three 
point contact plane. It divides the grain volume and the actual grain 
surface into two parts. Every three point contact plane has an extended 
three point contact plane counterpart. Since a grain may have vast 
numbers of contact points, especially when there is a matrix, many extended 


three point contact planes which can be formed may be coplanar (Figure 1). 
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Figure 1 Ellipsoidal Grain and Associated Contact Planes 


For a given surface grain, there is one particular extended 
three point contact plane (or group of coplanar. planes), the minimizing 
plane, which minimizes the volume of the part bearing all or the greatest 
portion of the apparent grain surface (1.3.0). This part bearing most of the 
apparent grain surface is termed the top part, the other, the bottom part. 
1.3.2 Three Contact Modes | 
The center of gravity of a grain is that point through which 
the resultant attraction of gravity, acts regardless of the grain's 
position. If the grain could be suspended or poised from this center, it 
would be in equilibrum in any orientation. Employing the grain assumptions 
in 1.3.0 the center of gravity is located in the center of each grain. 
With reference to 1.3.1 it may be seen that the center of gravity 
may lie in the top part, bottom part or on the minimizing plane. These 


represent the three modes of contact a surface grain may have (Figure 2). 
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Figure 2 Ellipsoidal Cross-Sections Illustrating 
Three Contact Modes 

The single point contact class may be considered a special 
case of contact mode 1, where the entire volume of the grain is the top 
part. For the two point contact class, it is conceivable that the grain 
is balanced and resting on two point contacts (e.g., a spherical grain 
resting in a V-shaped channel). 

It is likely that these two special mode 1 cases are rarely 
observed in nature for such reasons as grain surface roughness, instability 
and presence of matrix. 

1.3.3 Graphic Experiment 1] 

For an ellipsoidal grain the axial dimensions of its outline 
trace (1.3.0) and minimizing plane are dependent basically upon its k 
values, orientation and contact mode. Subsection 1.3.3 investigates 
this complex topic with a number of simplifying assumptions and a small 


graphic experiment. 


vi 


ia : 
| 
en an fala aa 
‘ ‘ ” 
/ ., f 
- — * 
4 ‘ F 
a 
en 
* : 
. Z 
ors ——— ah Sy 


f ; Vig )j d UY : “of PP | re 
atencom aa 
rin’? pad of sreaq > up Sen err ine 


iweheu Lit ota Aa ‘hy? aii ® 
wait) 08 7 7 ity nein 


red vf y ty itt freq stan ing progeat se 
nia, tails Ranga ART AMIS en 
ne: pay shines , ee) aaa, rae ony vg av) 


Vig a be 


Wisi re Coho, § org ids tees, Ae Pre, ks Bees ol na thy ae 
‘ 4 15 ae 7 7 


) TEGOIA SeMHRHLCN Ling & ney He’ zon hes 7 vaste Ty aes : 

5 Je i ‘ g Mi ' 

: ot o 
Tr 
ae paren ova stains ie IE 
- ; S« - io 
rH J AERTS : ian 1 ; nie , h Ke 7 ; 
ane sue - Mo fro) eneirs -Delas a Wie Pehiws rai a offi i. 
: et 


of, =* ys 4 zi 


ai 
a 


Laaieat det Yrikatang ni irbs dhs ace cL fix ue Th mee i? eh? 
taste dnb nega tig: ior jan ry raft we Adin y a 


$ ha) en i eS) a (7 a ae T “tay ‘g fsa 7 yeu - ia Lore 


hs ’ » 
om ' as Caen 
5 Ge) 
a ; , ay 
it We) 7G 
: ’ 
<< rar) ee ee : 
+ > 7 
> He: = ui} ral 


Three grain shapes are investigated, a sphere (k)=k,=1) and 
two ellipsoids (k)=k,=0.75 and k,=k,=0.55). The sphere is included 
because it is the simplest case; the ellipsoids, because their k values 
are the upper and lower limits of those commonly observed. 

In all cases the ellipsoids are oriented such that the A axis 
is perpendicular to the line of sight associated with the apparent 
surface (1.3.0). This axis functions as an axis of rotation for the B 
and C axes which it intersects in the center of the grain (center of 
gravity). The B axis forms angles (Or 622250 45°06 7250,900) with this 
line of sight. 

The spherical and oriented ellipsoidal grains are analyzed 
with respect to their behavior in each of the three contact modes. For 
contact modes 1 and 3 the center of gravity is arbitrarily situated 
0.24 X B above and below respectively, the minimizing plane. In all 
three contact modes the actual grain surface associated with the bottom 
part is never visible (buried) and the minimizing plane is normal to the 
line of sight. 

Combining these different grain k values, orientations and 
contact modes many variations are possible. These are examined by halving 
each grain so that the resulting planes are perpendicular to the axis of 


rotation. This sectioning procedure produces three distinct plane 


ni i 
minor axis _ 9 75 


minor axis _ 79 
7 major axis 


: ), and two ellipses ( 
major axis 


figures, one circle ( 
and 0.55). For the ellipses, the major axis corresponds to B, the minor 
axis to C, of the parent ellipsoid. Each of these figures have two 
measurable dimensions whose values depend upon the specific characteristics 


of the parent grain: 


1) The minimizing plane dimension (MPD) is the line resulting 
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from the intersection of the minimizing plane and the plane 
produced by sectioning (Figure 3). 

2) The apparent surface dimension (ASD) corresponds to the 
minor apparent axis b, of the grain's outline trace on the 
apparent surface. The ASD is determined with reference to 
the line of sight which passes through the center of gravity. 
The length of the ASD is the sum of the length of those two 
lines on the plane that are perpendicular to and lie on 
either side of the line of sight, are in the top part or 


on the MPD, and have maximum possible length (Figure 3). 


line of sight Figure 3 The Minimizing 
Plane Dimension (MPD) and 
Apparent Surface Dimension 
(ASD) of Oriented Ellipse 
Minor axis _ 0.75 


Bee grain - Major axis 


major axis = B axis of 
parent ellipsoid 
minor axis = C axis of 
parent ellipsoid 


c.g. center of gravity 

y angle in degrees 
between the major axis 
and the line of sight 
passing through the 
center of gravity. 


“ bottom part 


In the small graphic experiment (Figure 4) the major axes of 
the circular and elliptical sectioning planes are made equal. Because 
the axes of a circle are equal any change in their orientation makes 
no dimensional difference. For each case the rp percent is calculated 
(Table 1). Since a surface grain's apparent grain surface and 

MPD 


minimizing plane are both elliptical figures, in each case the ASD 


percent quantitatively relates b to the corresponding axis of its 


re 74 
te sun 
> 
oor?) 40 aL ie tie He ; ie] riGeel'a mai 
ay i259 Bntl yt ie) tu. a elem ae 
wh ee arey “a a) A ay? 5 tle an 
po te. uit peal ee ie 
ant 
a? “tor iJ i uy &F oa paz "4 wlan ou 
| + : Ahk | : 
wm oF f<fe hleqwe eis ah er ats 2/1} rr) ve ne i soy wd 
| pret ae 
) 8a ; ; gia 4. Ghat gay 3p ad et aP ce 
Lc owt trian 209 WeteKGe CPA LA ty) Settd 
: cry - 
at 62 ee ie ee ~ ris Jit ta, Te) oat f. 
UTE ° i a } Ber Sw Bey 1 ret 4 y 
nie gare * oh : a ! 1 a) in 
sous TES bua! : : ) 
; EN GAT Oo eenk De it, 29 
ee Ba 
eae 0 - ACAS UNL AP 
’ biaraqet 420° 
y i 
H Ws 4 yan 
‘ rhil« ’ 7h 
ee WS 2 SARS 
iba es at GW 


bint ae wae Wd ede! a4 


ae 


tay Bq, 2 
aad Miuety wrt 


Yo avtep, solan 


aeines «| 


,_ Saw! aby Nae v~ 


bei. vi, 


ref 


a 
i | 


se 
satu Pe eh Kh, > 


‘i 


00 a Cie Snails 


a Vrit 


PNP CS a a yi 


betike @qtt nd 


54 4a | a Bn? nes ore 


yes wid | 
zat} 


seth Fawshe nih pi 


vy 


si 
ore 


4 


WM: SS08RF VA Cey ys se eae py, 


er 
4 of 0 ie pas Ta 
ae 


: ‘Vadeit tL (jj BOAR. 7 


eee at c* 


ya] i nae he § 


10 


Saunhl4 auerg LedOLydti Ly pue uepnoulty 40 
JOLAGYagG SpOW Oe 2U0D *{ JUeWLUadxXy OLYdeug py SsunbLy 


€ 9pow 
39e3U0) 


asv°ddW 


2 2PpowW 
392 U0) 


dW 


L apoy| 
49e U0) 


eiainioia eietelwicay etn aes 


(SS°0) asditta (§Z°0) asdtt 13 BLIAL 


1] 


lL. & ve 


i “Gr Gz 


G'h2 
62 
62 
o¢ 

Ge Sc 


cg 
Lv 
G°6E 
cf 
Sg: 6c 


a 
60 
GeV 
€€ 
G°6¢ 


G°9C 
State 
Sve 
Gaac 

ve 


cS 
6P 
St 


GOP 


G°6€ 


ev 
S° LP 
Se asie 
ge 
S55 


S (9b 
Goce 
S295 
G Se 

ve 


os 
S°6t 
9t 
S°Ot 
G6£ 


oS 
GOS 
SOV, 
NY 
0°6€ 


OOL 


OOL 


9°48 


Saco 


9P 


G*2S 


J4ed woz Og 
UL APLARUY 
$O uazUue) 

€ 


aueld 

BuLZLwLuLy 

uo AZLAeUYN 

$O uaqua) 
c 


queg doy 

ut AP LARUYD 

$O UazuUay 
L 


(uu) 
adW asy % 


GG°0) asdtLL3a 


( ww) ( Wu!) 
ddW Gsv % 
GL°0) asdtiL3a 


GS¥ ‘ad Ss sys (Saeu6aq) apow 
GdW BL DALY A 39e3uU0)9 


| qUawLuadxy oLydeug Jo sziNnsay 
L 9Lgel 


12 


minimizing plane. 


For contact mode 1 grains, the mean aa percent decreases with 


minor axis 


TOMEI TIE value (Table 1). For the ellipse (0.55) in 


decreasing 
particular, minimum percent values are achieved when the major axis 
parallels the MPD. Brief consideration of Figure 4 suggests that as the 
center of gravity - minimizing plane distance increases so the Tey Percent 
values tend to decrease. 

For contact mode 2 grains, the mean Igy percent is nearly 100 
for the three planes, in all tested orientations (Table 1). The rey per- 


cent for the circle is 100 but decreases slightly with decreasing 


minor axis 


: MPD 
major axis value (ellipses 0.75 and 0.55 have mean ASD percent values of 


99.36 and 96.76 respectively). The minimum eae percent is reached for 
the ellipse (0.55) at an orientation of 45°, 

For contact mode 3 grains, the Fen Percent is uniformly 100 
(Table 1). Careful examination of Figure 4 diagrams reveal that if the 
center of gravity is closer to the minimizing plane, the Fey percent 
value drops slightly for the two ellipses (particularly ellipse 0.55), 
but not for the circle. For ellipses 0.75 and 0.55 noticeable reductions 
in rep percent begin when the center of gravity is 0.10 and 0.17 X B 
respectively, below the MPD. 

Figure 5 displays Table 1 data as well as a few additional 
contact mode 3 data points. Table 2 summarizes the trends in Table 1 and 
Figure 5. Table 2 is based entirely on results from ellipses whose major 
and minor axes correspond to parent grain axes B and C. These results are 
equally valid if the rotational axis perpendicular to the line of sight 
is B or C, rather than A. While this seems to extend the range of 


application for the results, the initial experimental assumptions prevent the 
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determination of the major and minor axes for both the minimizing plane 
and trace outline, except for a few special situations. The hep Percent 
values (Figure 5), however, and observed trends (Table 2) not only compare 
axial lengths but, indirectly, shapes and areas of the minimizing plane 
and the apparent grain surface which is bounded by the trace outline. 


This is particularly evident for the spherical grain. When in contact 


MPD 
ASD 


plane and the apparent grain surface are circular the diameter of the 


mode 1 the percent is 87.6 (Table 1). Since both the minimizing 
minimizing plane is 87.6 percent of the apparent grain surface. Thus 
the minimizing plane area is also much smaller. In contact modes 2 and 
3 the rep percent is 100 implying that for a spherical grain the minimizing 
plane and apparent grain surface coincide in shape and area in contact 
modes 2 and 3. 

Similar reasoning can be applied to ellipsoidal grains. Table 
2 and Figure 5 reveal that as ellipses decrease in aehvaneesinne the 
greater the potential for orientation to have a minimizing effect on the 
Tey percent. Figure 5 demonstrates that even for an ellipse (0.55) in 
any orientation, the minimizing effect is less than 9 percent for mode 2 
(re percent for a circle is 100 and for an ellipse (0.55) at 45° is 91.1) 


and quickly vanishes in mode 3. In contact mode 1, the effect of small 


minor axis 


; ; Hanes 6, MPD ; 
: : mem Percent 1S 
major axis values and orientation on minimizing the Asp Pe 


greatest. 


The above discussion provides evidence that ellipses, including 


minor axis 


—————— value is equivalent to common k values, are 
major axis 


those whose 
Sixt : ’ . MPD : 
similar to circles with respect to their Asp percent behavior for a 


given contact mode. In contact mode 1, for both spheres and ellipsoids 
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the apparent grain surface area and axial dimensions are greater than 
those of the minimizing plane. For contact mode 2, the apparent grain 
Surface - minimizing plane correspondence is exact for spheres and very 
close or exact for ellipsoids, especially with high k values. For contact 
mode 3, coincidence is always exact for spheres and eventually for any 
ellipsoid (the greater the center of gravity - minimizing plane distance, 
the more probably coincidence is exact for any k value and orientation). 


1.3.4 Graphic Experiment 2 and Description of Overlapping 
Grain Surface (0GS) and Embedded Grain Surface (EGS) 


Where a material (1.3.0) has an approximately planar surface 
consisting of an assemblage of grains in one or more contact modes, this 
material may have two distinct types of surfaces; overlapping grain 
surface (OGS) and embedded grain surface (EGS). These are defined in 
terms of which contaerinndes are associated with the surface. Graphic 
experiment 2 reveals their characteristics for a material consisting of 
ellipsoidal grains. 

An OGS consists of an assemblage of grains belonging to all 
three contact modes. The characteristics of this surface may be 
appreciated by first considering all grains to be spheres. Following 
this discussion the more complex ellipsoidal grain situation is examined. 

It is assumed that a material consists of identical spherical 
grains and the minimizing plane of each surface grain is coplanar with 
the surface plane. Since this is an OGS, there are surface grains 
associated with all three contact modes. 


Assuming that there are equal numbers of grains present at 
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any center of gravity - surface plane distance (A) up to and including 
One radius, then all spheres within this A range touch the surface plane 
and are represented by their respective minimizing planes. Beyond one 
radius, the spheres abruptly cease to touch the eases plane and are 

no longer represented on it. 

Another consequence of the equal numbers assumption is that 
contact modes 1 and 3 grains are equally numerous. Contact mode 2 grains 
are very rare in comparison because they represent the special case 
where A is zero. 

In the OGS spherical grain case, orientation effects are non- 
existent. For ellipsoids this is not so. Since the ellipsoid may have 
three different axial values, grain contact with the surface plane is 
dependent on both A and grain orientation. In the following discussion, 
assumptions similar to those in OGS spherical grains are employed, the 
primary exception being that the material consists of equi-sized 
ellipsoidal grains. 

Graphic experiment 2 demonstrates the dependence of surface 
plane representation of ellipsoids on A and orientation, In this 
experiment, which is analogous to the one in 1.3.3, two ellipsoids 
=k,=0.75 and k 


ite? ] 
that the A axis is perpendicular to the line of sight associated with 


(k =k=0.55) are investigated. They are oriented such 
the apparent surface. This axis functions as an axis of rotation for 
the B and C axes which it intersects in the center of the grain (center 
of gravity). The angle y between the line of Sight and major axis has 
five values (0052250 ,45°,6/ 250590) These oriented ellipsoidal 
grains are examined with respect to their behavior in contact mode 3. 
The grains are situated so that for a particular orientation only one 


point on their surface touches the surface plane. 


grithetpat tea ay on (A) ee suivant 

unig 40d ed ; oa & bi) sehen 
Ys, 

ono DHoPSA: <9 Pe payne eve saya Niue 

awh toin. tee A eae 4 Ou ain au sea 


~ af m8 


ty 


: Siu -yadimin “ley ana iq auneoehgl 


“M 


ghinag © ube. Snes. suivante 9 Satnotiain 


sah the ies nit: tude) bee es G oa 
2s Rates ete cet q is apart viet t yn i ” 
¥ F rn ; i i ta 


+ < oe 


wi" iS, s at Se LS KC { le nST } oa ae Pa As Giga 1ni i i), 


A yee bhocarTits 447 e307 a a, “ue mys 7 


27 STR, s wt ie OAT ATP ar tnd ates iy ear ban (4 


ra = 


yaad. pawl To} mets ol mipyhie: iio heb ; 
¥) .tagbtowe ath 2ahee Leena maa al seni oF = 
_ baste vey Ww tee tenes Leinspie oi ‘tent be 


cath 
ue he ‘4 ¢; La) ? 
or rusk 


isemiyn Vo samibngest odf 2988" rede =) ee ret oe 
Tinh at eee ae Ie Ae. seen ss ty led 


vf 


oh hy 
Cw G?atiio wet ,f:2.1 at #50, o0) wo! wil br oF fo) teats 


ye pak feriteg vtnresit . badnghtene int eRe 7) | brie ahi ce s 
ATP BOntocKes shal «Sy oF ating Ye fu tinidery ener tai batt a an 
: e ; a 3 vd he 


- . ve = } 
“a ain nghete: te. 27 ae ne Be sO (SSA CUE ai vil ais ‘inal ‘wid Ki a ny 


t j ‘ at | 

aie hin Pree yt? yatnns bis FERRARA Tt i sind } ye a La. Natt i 
4 ; : J s : ey _ 
7 ee pine Bar \' i MT Op ori Ts Les pe! * en bnrenig Wa ¢ be Pe 


7 


Gabi ‘eat al Dn) dere pen deat 8,8) Re iy) tu ey RTT, i j 


a a 7 pk 7 ; ; - = 
Fads O 1) “Oe RIE, ATSAG ee kt AI gat trae ave pniesy Hert 
gs 


MET” Nes Oe" hy ‘neddig vans asteny wt 


~ 


ees Spp3 = sata fr aap tan : 


18 


As in 1.3.3 each ellipsoidal grain is halved forming two 
ver . MmipOrraxis: = C. 2 
distinct ellipses Gakonsamis > B= 0.75 and 0.55). Each type of 
ellipse is drawn (Figure 6) in its five different orientations and with 
only one point on its boundary touching the surface plane. A is then 


measured for each of these orientations (Table 3) and plotted (Figure 7). 


Table 3 Results Graphic Experiment 2 


y (degrees) Ellipse (0.75) Ellipse (0.55) 
. A=%B A=%B 
90.0 S75 28.0 
67-5 40.0 310 
45.0 43.8 40.5 
22-5 48.6 47.6 
0 0 50.0 


Figure 7 shows that as the orientation of both elliptical 
types changes from y=90° (a=0.5 x C axis) to y=0° (a=0.5 x B axis) a 
increases. Closer scrutiny of the two curves in Figure 7 uncover their 
nature as the limits of total surface plane representation for a given 
elliptical shape. merninees with A and y values less than or equal 
to their respective curves, touch the surface plane and therefore are 
represented on it. For combinations exceeding curve values the converse 
is true. 

If it is assumed that there are equal numbers of ellipses at 
any.given y, and all orientations of B about A are equi-probable, then 
Figure 7 curves correspond to surface plane representation frequency 
curves, the frequency of. representation being dependent upon A andy. 
Figure 7 illustrates that representation on the surface plane is total 
up to a critical 4 (50 percent C), then steadily declines to a point of 


zero representation (50 percent B). Furthermore, the more circular in 
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Limits of Total Representation of Elliptical Grains 
on a Surface Plane (refer Table 3) 
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form the ellipse the greater the A before orientation - induced decline 
is initiated. As in 1.3.3 identical findings are produced if the axis 
of rotation is B or C instead of A. 

The results of this graphic experiment also hold true for the 
parent ellipsoids. More generally, it is expected that similar A-orien- 
tation dependent behavior is exhibited by ellipsoids in which any orien- 
tation is equi-probable. If so, the more spherical the ellipsoidal form, 
the greater the A (maximum 50 percent A axis) before orientation - surface 
plane representation effects are initiated. 

The above observations are based on an experiment using only 
contact mode 3 ellipsoids. It is clear that for both contact modes 1 and 
3 grain representation on the surface plane is dependent upon identical 
A-orientation combinations. Certain grains do not touch the surface plane 
because they are buried; identical potential contact mode 1 grains are 
not possible because they are not supported, assuming that the minimizing 
plane of all grains of this surface are coplanar with the surface plane. 
Therefore, the surface plane representation frequency curves for contact 
modes 1 and 3 are symmetric about the surface plane. 

These statements remain applicable if a material bearing an 
OGS contains numerous various Sea SAAS, Instead of the surface 
plane being touched by identically shaped grains in different a-orientation 
combinations, the surface plane is shared by representatives of each size. 
This complex situation closely corresponds to OGS observed in nature. A 
familiar example is the surface of exposed gravel bars which usually 
consists of ellipsoidal shaped pebbles and cobbles. Contact mode | clasts 


rest upon and thus overlap contact modes 2 and 3 material (Figure 8). 
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Figure 8 An Apparent OGS 
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An EGS consists of an assemblage of surface grains belonging 
to contact modes 2 and 3. Assuming that the minimizing plane of each 
Surface grain is coplanar with the surface plane and the material consists 
of numerous size fractions of ellipsoidal grains, She behavior of the 
contact mode 2 and 3 surface grains is identical to that for OGS. Thus 
contact mode 2 grains are rare in comparison to contact mode 3 grains and 
the surface plane representation of these latter grains is dependent upon 
A and orientation. 

Both artificial and natural materials may possess an EGS. 
Usually these materials have a matrix. The surface of the matrix is 
generally equatable with the surface plane of the EGS. 

Artificial materials possessing an EGS, typically are indurated. 
A familiar example of EGS is a paved road surface of sand and gravel 
mixed with asphalt. During paving the asphalt is soft and the grains 
are embedded by rolling. Subsequent to hardening some grains may become 
contact mode 1 grains due to the disappearance of some of the surrounding 
asphalt matrix. This creates an ephemeral OGS, since these grains are 
easily dislodged. 

In nature, gravels which have fine sand, silt or clay matrixes 
can maintain near vertical faces, as is often the case with terrace 
gravels. Investigation of these faces reveals that they are EGS. The 
explanation for this is quite similar to that for road surfaces. If a 
grain on a vertical face loses some of its supporting point contacts 
(matrix or similar sized grains) and becomes a contact mode 1 grain, the 


grain weight is unsupported and it tumbles out of its mold (Figure 9). 
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Figure 9 An Apparent EGS 
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1.3.5 Comparison of Apparent OGS, Apparent EGS and 
Thin Section 

If a randomly chosen plane is extended through a material 
(1.3.0) consisting of ellipsoidal grains, the frequency of grain repre- 
sentation on this plane is the same as discussed in 1.3.4 (O0GS Ellipsoidal 
Grains). If the material is cut along this plane, a planar surface 
results consisting of intersected grains and matrix, if present. Referring 
to 1.3.4, the cutting plane is equivalent to the surface plane, and the 
intersected surfaces of the ellipsoids are equivalent to minimizing planes 
coplanar with the surface plane. This planar surface can be considered 
a thin section. 

For OGS, the surface plane contains minimizing planes of all 
three contact mode grains. As discussed in 1.3.3, the minimizing plane 
and apparent grain surface of contact mode 1 grains do not coincide, 
Coincidence is far more probable for contact mode 2 and particularly for 
contact mode 3. Since nearly half of OGS grains are in contact mode 1 
(1.3.4), and they lack coincidence, an apparent OGS is not equivalent to 
a thin section corresponding to the surface plane. More generally, since 
all thin sections cut from this material should have equivalent grain 
representation, the apparent OGS is not equivalent to any thin section 
surface. 

For EGS the surface plane contains minimizing planes of only 
contact modes 2 and 3 grains. Coincidence between the apparent grain 
surfaces and their minimizing planes is very good, particularly as A 
increases (1.3.3; 1.3.4). This seems to ensure apparent EGS comparability 
with its associated thin section surface, though it may be argued that 


the center of gravity of grains represented on a thin section are on 
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either side, or on the cutting plane. In contrast, EGS grains have their 
center of gravity on or to only one side of the surface plane. The 
effect of this is negligible for the following reason. As discussed in 
1.3.4, surface plane representation frequency curves for an OGS are 
symmetric about the surface plane. This also applies to thin sections. 
Since contact mode 1 grains are non-existent in an EGS, the 'gaps' are 
filled by contact mode 3 grains. On the basis of surface plane repre- 
sentation symmetry, it is probable that the replacement contact mode 3 
grains will pemreneccentc’ on the surface plane in a very similar 
manner. Therefore, the net result is an apparent EGS which is equivalent 
to any thin section cut through the material. 


1.3.6 Grain-Size Measurement on Apparent OGS and EGS 
and Thin Section 


If the apparent surface of a material (1.3.0) is used for 
grain-size analysis it is necessary to assign linear dimensions to the 
trace outlines of the surface grains. There are a variety of measures 
which may be employed (Kellerhals et al 1975). 

Apparent Axis Measurements 

a) major apparent axis a 

b) minor apparent axis b 

c) some combination of a and b 

Chord Length Measurements 


a) chord length of a grain falling along a predetermined line 
b) maximum chord length in a predetermined direction 


For an apparent OGS approximately one half of the surface 
grains are contact mode 1. The fabric may be anisotropic such that true 
axial values may be obtained directly by measurement of the apparent 
axes of contact mode 1 grains. At least one true axis must be parallel 


with the surface plane. 
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a) If only one true axis of the surface grain is parallel 
to the surface plane and its direction on the plane is 
fixed, measurement of the corresponding apparent axis, 
provides the true axial value. 
b) If one true axis of the surface grain is normal to the 
surface plane, then the other two are parallel to the 
plane. Association of the two apparent axial values 
with the corresponding true ones may be done on the 
basis of relative size, preferential direction or both. 
Particularly for contact mode 3 grains it is improbable that 
an apparent axis corresponds with a true one. Because these grains are 
the primary constituents of an apparent EGS this fact must be recognized 
by all grain-size analysis procedures utilizing this apparent surface. 
Thin section surfaces and the apparent EGS of a material are essentially 
equivalent (1.3.5), and it is reasonable to assume that thin section 
grain-size analysis may be applied to an apparent EGS. 
1.4 Sampling Procedures 

1.4.0 Introduction 

Section 1.4 examines all possible sampling procedures which 
can be used on thin section surfaces and the apparent EGS of terrace 
gravels. When the constraint imposed by the necessity for geometric 
equivalence of a sampling technique to bulk sieve analysis is taken into 
account, very few of these sampling procedures remain applicable. 

Kellerhals et al (1975, p.80) state that all grain-size 

analysis procedures in geology may be classified according to three basic 
steps and associated choices. These steps are: 

Step I Collection of sample 

Step II Assignment of linear dimensions to grains 

Step III Allocation of frequencies 


In the following discussion this clear and comprehensive system of steps 


and choices is adopted. 
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1:4.1. Step I Collection of Sample 


The basic choices associated with collecting a sample are 
fixed by the number of dimensions (i.e., zero, one, two, three) of the 
total sample volume (sum of the volumes of all grains in a sample) 
determined by the experimenter, and the number of dimensions dependent 
On grain size. 

1) Volumetric (Bulk) Sample: The sample consists of a 
volume of the material under investigation. Three 
dimensions are predetermined by the operator. 

Standard sieve analysis uses this method. 

2) Areal Sample: Every grain in a given area is examined. 
Two dimensions are predetermined by the operator. 

3) Transect Sample (line counting): A line is placed 
Paeoueh the material. The operator predetermines one 
dimension. 

4) Grid Sample (point counting): Only dimensionless 
points are predetermined within the material. 

Only methods 2), 3), and 4) are applicable to thin section analysis. 


1.4.2 Step II Grain-Size Measurement and Step III 
Frequencies . 


In 1.3.6 various methods are presented for assigning linear 
dimensions to trace outlines of surface grains on apparent surfaces. 
These methods are also applicable to trace outlines on thin sections. 
These grain-size measures, as well as other methods of assigning linear 
dimensions to grains, constitute Step II (see Kellerhals et al, 19/5, 
p.80). This step is Nees eary for computing frequencies, Step III, since 


frequency determination requires that each grain be a member of a size 


interval. 
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1) Frequency by Area: The frequency of each size interval 
is expressed as the percentage by area of the original 
sample falling into the interval. 

2) Frequency by Length: The frequency of each size interval 
is expressed as the percentage by length of the original 
sample falling into the interval. 

3) Frequency by Number: The frequency of each size interval 
is expressed as the percentage by number of the total 
number of grains in the original sample that fall into 
the interval (Kellerhals and Bray, 1971, p.1169). 

4) Frequency by Weight: The frequency of each size interval 
is expressed as the percentage by weight of the original 
sample falling into the interval (Kellerhals and Bray, 
1971, p.1169). 

Frequency methods 1), 2), and 3) are applicable to the analysis of 
apparent OGS and EGS and thin sections. 


1.4.3 The Concept of Geometric Equivalence and 


Its Importance 


A number of different sampling procedures arise because of the 
alternative methods which exist in Steps I and III. It is necessary to 
determine whether these sampling procedures are geometrically equivalent. 
"Equivalent sampling procedures, applied to a homogeneous and isotropic 
deposit result, on average, in identical size distributions" (Kellerhals 
and Bray, 1971, p.1166). 


There are several important issues associated with the concept 


of geometric equivalence. 


1) Although many sampling procedures have been used, questions 
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dealing with their geometric equivalence are rarely 
adequately explored. This has resulted in the comparison 
of non equivalent data and in other cases the omission 

of important results because differences in procedures 
could not be evaluated (Kellerhals and Bray, 1971, p.1166). 

2) Customary bulk sieve analysis (volume-by-weight sampling 
procedure) is central to the problem of equivalence 
because most grain-size analyses employ this procedure. 

As such, most accepted theories of sedimentation are 
based on data derived from bulk sieve analyses. Thus, 
it is critical that other sampling procedures, such as 
those applicable to thin section grain-size analysis, 
yield comparable grain-size distributions to those that 
would be achieved by this standard sampling technique. 

3) Recognition of the problems associated with geometric 
equivalence have lead to geometric formulations of this 
concept and its implications. Ensuing geometric arguments 
reveal that some sampling procedures are directly equiva- 
lent while others must undergo conversion to achieve 
comparability. 

1.4.4 Determination of Conversion Factors for Steps I and IIT 

The following discussion presents conversion factors for Steps 

I and III. Because of the importance of bulk sieve analysis, conversion 

factors are determined relative to volume (Step I) and weight (Step III). 
Since many basic assumptions used to derive these factors are essentially 
those found in Kellerhals and Bray (1971) the derived factors are 


identical to those used by them. 
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Assumptions 

A) A homogeneous, isotropic specimen material (1.3.0) 
consists of densely packed cubes of constant specific 
weight. These cubes may be referred to as grains 
(Figure 10). 

B) There are n different size fractions each consisting of 


One or more cubic grains of linear sizes D> apy i) 


a Sie ener s 


Dn 
C) The specimen volume V can be divided into constituent 
volumes Vi> Vos V35 Ase ee each representing the total 
volume of the cubes of the corresponding size fraction. 
These volumes form a ratio Vz:Vo:V3: tas Vee 
D) In a volumetric (bulk) sample all three dimensions are 
predetermined by the operator (Kellerhals et al 1975). 
Thus, none of the sample dimensions are dependent on 
grain size. If V is sampled volumetrically the ratio 
of the volumes associated with each size fraction is 
expected to be equal on average to the ratio V1 2V5:V3: ers 
Vie 
E) A cut parallel to a side of the cube produces a 
smooth planar specimen surface S, of area A, in which the 
exposed area of each surface grain is proportional to the 
square of its linear size 
Employing these assumptions six conversion factors are deter- 
mined. 
1) Areal Sample to Volumetric Sample 


It is assumed that the grains of each size fraction cover a 


total area (a) on S, such that these areas form a ratio A, iapidzi ... a 
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Total No. 


Total No. 


Linear 
Size D 


in Sample 
Surface 


in Sample 


Volume 


Weight 
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Grain 
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Figure 10 
Sample of Densely Packed Cubes of Three Sizes 
(from Kellerhals and Bray 1971, 
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which is equal to the ratio of the corresponding volumes V,2V5:¥3: yas 
_ (see Kellerhals and Bray, 1971, p.1175). 
Areal collection involves the measurement of every grain on 
S. Volume values may be derived from this a Sample if the volume of 
every surface grain is calculated. 
The volume of each size fraction is a, X Ds a5 X Dos az X D3 


ot Moe Do: But the volumes of size fractions collected areally will 


n 
not necessarily form a ratio equal to Vj 2V5:¥3: . me Equivalence 
is achieved by multiplying each size fraction by r. 

Lae a De 1 
Thus a, x D, x D, > Gy X dD, X D, a3 X D3 X D, Sse. A,X 0, x D, 


a5 geo eee an 


V1 iVo:V9: Eee YS 


2) - Transect Sample to Volumetric Sample 

It is assumed that for any transect which is placed on S, 
the ratio formed by the combined lengths of those grains in each size 
fraction which the transect touches, 1, >15513, ns 1 is Tyrlgtl3: 


and is equal on average to Vy: Vo:V3: oe | 


n 

Transect collection involves the measurement of every grain 
which it touches. Volume values may be derived from a transect sample 
if the volumes of collected grains are calculated. 


The volume of each size fraction is 1, x DS, 1, x tee x 1, x D>» 


‘ 1, X of But the volumes of size fractions collected by transect 
will not necessarily form a ratio equal to V,2Vo:V3: ae wes Equivalence 


is achieved by multiplying each size fraction by ne. 
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3) Grid Sample to Volumetric Sample 

It is assumed that for any origin point and orientation of a 
grid on S, the ratio formed by the number of those grains in each size 
fraction falling under grid intersections, 9129999 oe is 
G72Go:93i +++ GD, and is equal on average to VyiVoiV3... va 

Grid collection involves the measurement of each grain falling 


under each grid intersection. Volume values may be derived from a grid 


sample if the volumes of the collected grains are calculated. 


: on eae 3 3 3 
The volume of each size fraction is 9, x D> Jo Xx DD 93 X D2, 

- g, X D°. But the volumes of size fractions collected by grid will 
not necessarily form a ratio equal to Vi:Vo:V3: can “ee Equivalence is 


achieved by multiplying each size fraction by 5 
Sie a ls ete Sel Page Soha 
Thus g, x D, x De > Go X Do x 9a x D3 x D° So aes Gy X D., x De 
] 2 3 n 
=91299'93: «+. Gy 
ViiVo:V3: Bi vi 
In frequency by weight the frequency of each size interval is 
expressed as the percentage by weight of the original sample falling into 
the interval (1.4.2). If the grains have a constant specific weight, 
weight frequencies are identical to volume frequencies (Kellerhals et al, 
1975, p.80). Thus for size fractions in the specimen material: 
V1 iVo:Vo: cre ve = Ws 2Wo:W2: ‘2 W, where W is the weight of the volume 
of a given size fraction. 
4) Frequency by Area to Frequency by Weight 
To convert frequency by area to frequency by weight the calcu- 
lated area of each size fraction must be converted to a volume. This is 


performed by multiplying the area of each fraction by its corresponding 


linear size, D. The resulting volumes are proportional to weights. 
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5) Frequency by Length to Frequency by Weight 

To convert frequency by length to frequency by weight the 
observed grain numbers of each size fraction must be converted to a 
volume. This is performed by multiplying the length of each fraction 
by its corresponding area, De. The resulting volumes are proportional 
to weights. 

6) Frequency by Number to Frequency by Weight 

To convert frequency by number to frequency by weight the 
observed grain numbers of each size fraction must be converted to a 
volume. This is performed by multiplying the number of each fraction 
by its corresponding volume, D°. The resulting volumes are proportional 


to weights. 
1.4.5 Combined Conversion Factors for Nine Sampling Procedures 


Nine sampling procedures emerge when methods applicable to thin 
section analysis, within Steps I and III, are combined. For each 
procedure, the conversion factors derived for densely packed cubes in 
random arrangement (1.4.4) may also be combined. These combined conver- 
sion factors e.g., 1, D, u can be used to adjust the value of each size 
fraction as determined by any of the nine sampling procedures, to its 
volume-by-weight equivalent. More generally, the results of any procedure 
can be transformed so as to be equivalent to any other (Table 4). 

Three sampling procedures, area-by-area, transect-by-length and 
grid-by-number have combined conversion factors of one with respect to 
volume-by-weight, and thus are geometrically equivalent to this standard 
and to each other. Identical conclusions were reached by Kellerhals 


et al (1975). 
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Table 4 
Thin Section Sampling Procedures 
and Combined Conversion Factors 


Combined Conversion Factor: 


Steps and Sampling Panuersion. fo 
Methods Procedure Volume-by-Weight 
2a LETC) area-by-area tx D=1 
] al 
1(3) = ITI(1) transect-by-area “5 X D FF 
D 
. ] es 
14) = 111(1) grid-by-area —, xD=z- 
3 2 
D D 
1(2) - I11(2) area-by-length aX pate B 
Ti3)ie Fi ¢2) transect-by-length 5 X Ng = ] 
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1.4.6 Geometric Equivalence and Equivalent Grain-Size Measures 


Most investigators agree that for commonly occurring grain 
shapes, the square mesh sieve size Ds and the intermediate axis B (1.2.4) 
are both acceptable and almost identical measures of actual grain size 
(Leopold, 1970; Kellerhals and Bray, 1971). Thus, for the three sampling 
procedures in 1.4.5 which are geometrically equivalent to bulk sieve 
analysis, measurement of the B axis of the collected grains should yield 
a grain-size distribution closely similar to the one produced by sieving 
the specimen. 

Under certain conditions, there may be direct axial correspond- 
ence between an apparent and true axis for grains on an OGS (1,3.6). 

If the true axis is B and the corresponding apparent axis is measured, 
it seems likely that any of the three geometrically equivalent sampling 
procedures (1.4.5) may be used to derive a grain-size distribution which 
is nearly identical to the one achieved by sieving. 

For thin section and terrace gravel surfaces it is probable 
that the apparent axes of sampled grains do not correspond with true 
axial values (1.3.6). While sampling procedures yeometrically equivalent 
to sieve analysis can be employed on these surfaces, some technique must 
be devised to obtain "the feeri bution of actual arain size (A,B,C,D, ...) 
from the corresponding distribution of observed sizes" (an Deere) 
(Kellerhals et al 1975, p.82). These “observed sizes" refer to the linear 
dimensions assigned to trace outlines of grains on either an apparent 
surface or a thin section. 

1.5 Grid-by-Number Sampling of the Surface Layer of Exposed Gravel Bars 

1.5.0 Introduction 


Three surface-oriented sampling procedures, area-by-area, 
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transect-by-length and grid-by-number were shown to be geometrically 
equivalent to volume-by-weight (1.4.5). It is necessary to determine 
whether this theorized equivalence can be confirmed by sampling actual 
deposits. 

Much work of this kind has been performed using the surface 
layer of exposed gravel bars. In these studies of this OGS (1.3.4) 
most research has compared the grain-size distributions arising from 
grid-by-number and bulk sieve analysis. The problems associated with 
collecting a grid-by-number sample of this surface layer (especially when 
only the apparent surface is examined) and the equivalence of this 
technique to bulk sieve analysis have many implications to similar 
terrace gravel grain-size investigations. 


1.5.1 Comparison of Bulk Sieve and Grid-by-Number (B,) 
Distributions 


Kellerhals and Bray (1971) list four grid sampling methods 

which may be used on the surfaces of exposed gravel bars. 

1) A grid is established over the gravel surface and the 
grains immediately beneath the grid points constitute 
the sample (see Thornes and Hewitt, 1967; Kellerhals 
and Bray, 1971; McGinn, 1971). 

2) A survey tape is stretched across the area to be 
sampled using a set of regularly spaced points (e.g., 
footmarks) as grid points (Wolman, 1954). 

3) A grid system is formed by traversing the sample 
area collecting the stone immediately beneath the 
toe after one a more steps (Wolman, 1954). 


4) The operator walks along several parallel lines 
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Stopping after each stride, averts his gaze and reaches 
down over his toe with a finger and the first rock 
touched is picked up for measurement (Leopold, 1970). 
Once a grain was selected by one of the above grid collection 
techniques, usually one or more of the true axes were measured in the 
field (henceforth: subscript t signifies that the term's value depends 


upon the measurement of surface grain true or apparent axes; A By and 


‘Ee 
C. are terms denoting the field measurement of a grain's true A, B, and 


C axes; AB,» ABC, . 


derived from a grain's Ay» B> and C,3 ABC, is the triaxial mean of a grain). 


BC, are terms of denoting different arithmetic means 


Wolman (1954), and Kellerhals and Bray (1971) took B, as being comparable to 
the square mesh sieve size Des and used 50 to 100 axial values to compare 
bulk sieve and grid-by-number grain size distributions. McGinn (1971), 

as well as testing By. examined the consequences of using grain-size 
measures AB, > and ABC... 

Generally, the results from these comparisons indicated that 
grid-by-number, using By. or ABC, as a measure of grain size, yields 
grain-size distributions coarser than those generated by bulk sieve 
analysis procedures (henceforth: Arsg> AB 50° Bi sg: ABC, £9: BC sg? and 
Ci59: are terms denoting various medians of grid-by-number distributions 
based on field measurements of true axes; Avs AB, > Bis BC BC. and C, 
are terms denoting various means* of grid-by-number distributions based 
on field measurements of true axes; subscript sieve signifies that the 
term's value depends upon square mesh grain-size measurement; Di eves0 
is the term denoting the median of a grain-size distribution produced by 
sieving; D is the term denoting the mean of a grain-size distribu- 


sieve 
tion produced by sieving). While Wolman (1954) found that Bisg was 


*means used in this thesis are Folk and Ward unless otherwise specified 
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substantially greater than D Kellerhals and Bray (1971) concluded 


sieve50’ 
on the basis of 15 samples, B50 was only slightly greater than D 


sieve50 
(grains > 8mm). McGinn (1971, Table 4) found that grid-by-number (B,) 

and bulk sieve distributions (grains > 8mm) were significantly close in only 
43.1 percent of the 30 test samples. McGinn also utilized the Wilcoxon 

Test (Wilcoxon Matched Pairs Signed Ranks Test) to compare the statistical 
parameters (median, mean, standard deviation, skewness and kurtosis) of 

each grain-size distribution. The Wilcoxon Test indicated that grid-by- 


number (B,) and bulk sieve distributions were not statistically similar 


) 
except for skewness. As well, the Wilcoxon Test (one-tailed) revealed 
that B. > Derava: 

Most researchers agree that there is one major problem with 
using the surface layer of exposed gravel bars to test the equivalence 
between grid-by-number and bulk sieve sampling procedures. As Kellerhals 
and Bray (1971, p.1166) explain: 

"At low to intermediate stages virtually all sand and 

sometimes the finer gravel fractions are removed from 

the bed surface of a gravel-bed river. This results in 

a distinct pavement of the bed with a gravel layer one 

grain thick ... gravel beds commonly consist of two 

separate populations, the surface layer and the under-~ 

lying deposit." 

Kellerhals and Bray argued that while the surface layer is very important 
to hydraulic friction or initiation of bed movement studies, it theoreti- 
cally cannot be sampled volumetrically. This is because the gravel layer 
is only one grain thick, thus the thickness dimension of an intended 
volumetric sample (all three dimensions are predetermined by the operator) 
cannot be predetermined. This argument justifies the use of a surface 


oriented technique like grid-by-number but suggests it may be difficult 


to confirm the equivalence between this technique and customary bulk 
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Sieve analysis. If population differences are ignored completely the 
grid-by-number (B,) distribution is substantially coarser than that 
produced by sieving (Wolman, 1954). Both Kellerhals and Bray (1971) and 
McGinn (1971) eliminated material less than 8 mm Frei their bulk sieve 
analyses in what the former said was, " ... an attempt to compensate for 
the population differences between the surface layer and the underlying 
volume.". However, as discussed, both still found grid-by-number (B,) 
distributions to be coarser than those of bulk sieve. | 


1.5.2 Comparison of Grid-by-Number (B 


and (b,) 
Distributions 


) 
In some of the above studies, the apparent surface of exposed 
gravel bars (an example of an apparent OGS, refer 1.3.4) was sampled with 
grid-by-number. 
The surface layer of the exposed gravel bars was photographed 
so that the line of sight was approximately normal to the surface plane. 
A grid was superimposed either in the field or on the photograph (slide 
or print). In either case, the photograph always contained some means 
(rulers at right angles or the grid) whereby the surface plane could be 
more precisely normalized. With this done and the grid in place the 
trace outlines of the surface grains falling beneath the grid intersection 
points were measured. Generally, only those grains which were > 8mm 
were used (Thornes and Hewitt, 1967; Kellerhals and Bray, 1971; McGinn, 
1971). 
The problem remains as to how the trace outlines of the grid 
selected surface grains should be measured. Pashinsky (1964, p.279) and 
Leopold (1970, p.1358) observed that the C axis of the surface layer 


gravels were normal to the surface plane. Although this is an approximation 


in cases where there is imbrication (Johanson, 1963), basically the surface 
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layer can be treated as an example of a situation discussed in 1.3.6, 
case b) and 1.4.6, namely, there is a direct axial correspondence between 
major and minor apparent axes a and b and true axial values A and B 
respectively, for contact mode 1 grains of the apparent OGS. Generally, 
contact mode 1 grains can be recognized on the apparent surface because 
they tend not to be overlapped. Therefore, it can be expected that 
measurements of the b axis of the gravel bar contact mode 1 surface grains 
selected by a grid, would produce a grid-by-number distribution very 
similar to that of grid-by-number (B,) (henceforth: a, and b, are terms 
denoting the major and minor apparent axes respectively, of a surface 
grain associated with either an apparent or thin section surface; ab, 

is the term denoting the arithmetic mean derived from ay and b,3 A459? 
abi cg: bis are terms denoting various medians of grid-by-number 
distributions based on measurements of surface grain apparent axes; ays 
ab,» b, are terms denoting various means of grid-by-number distributions 
based on measurements of surface grain apparent axes). 

Kellerhals and Bray (1971) and McGinn (1971) compared the grid-— 
by-number distributions based on the measurement of B and by. The 
former advanced the tentative correction formula, bi 5970-888 i 59 based 
on 14 sample points. McGinn's results revealed that for 30 samples 
compared in this manner, the grain-size distributions were statistically 
equivalent in 93.3 percent of the cases. The Wilcoxon Test (McGinn, 
Table 5) disclosed that grid-by-number (B,) and grid-by-number (b,) were 
significantly close for median, skewness and kurtosis parameters. 

Kellerhals and Bray recognized that there was some disagreement 
between the two sampling procedures, but concluded that the results were 
closely equivalent. McGinn stated that these grid-by-number procedures 


produce significantly close results for median and mean values. 
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1.5.3 Comparison of Bulk Sieve and Grid-by-Number (b,) 
Distributions eres 


Both Kellerhals and Bray (1971) and McGinn (1971) examined the 
relationship between grid-by-number (b,) and bulk sieve grain-size 
distributions. Based on 1] sample points, Kellerhals and Bray advance 
Deca (grains greater than 8mm) = 1.0bi 55 + 5mm as a tentative 
correction formula. McGinn found that these two grain-size distributions 
were equivalent in 66.7 percent of his 30 samples. The Wilcoxon Test 
(two-tailed) indicated that the two distributions were significantly 
Similar only for standard deviation and skewness parameters and the 
Wilcoxon Test (one-tailed) confirmed that b, was greater than Lee ies 
The above results seem to conflict, since Kellerhals and Bray 


found D was greater than bi59 and McGinn showed that b, was 


ic 
Examination of 1.5.1 and 1.5.2 reveals that these 


sieve50 


greater than De es 


results are consistent with their respective studies. Since Kellerhals 


and Bray found Bi50 slightly greater than D UeateorleseeanG By 


50 
greater than bis0 in 1.5.2, it is quite possible that one should find 


sieve50 


Daan greater than bis0 (1.5.3). McGinn observed B. was greater 


than D in 1.5.1 but he found in contrast to Kellerhals and Bray 


sieve 


Bi 5g was equivalent to b+50 (1.5.2). Therefore, McGinn's finding that 


b, was greater than D 


; (1.5.3), is reasonable. Perhaps the C axis 
sieve 


of the surface layer gravels McGinn examined were more strongly oriented 


normal to the surface plane. 


1.5.4 Relative Coarseness of Grid-by-Number (B,) Distributions 
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The gravel bar experiments reviewed in this section indicate 
that grid-by-number (B,) distributions tend to be coarser than their bulk 


sieve counterparts. This may arise for two reasons: 
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Geometric Non-equivalence 

As discussed in 1.5.1 gravel beds often possess two 
separate populations due to surface paving, a surface 
pavement one grain thick and a finer underlying deposit. 
Given these conditions, the geometric equivalence between 
bulk sieve analysis and grid-by-number would breakdown, 
Since the latter employs exclusively the coarser surface 
layer. Although these population differences have been 
recognized and compensated for by only sieving material 
greater than 8 mm, this compensation may not be 
sufficient to overcome population differences. 
Non-equivalence of Grain-size Measures 

Discrepancies may arise due to the assignment of linear 
dimensions to grains. It has generally been assumed that 
the B axis is comparable to Des as a measure of grain-size. 
If population differences have been compensated for 
adequately, then the relative coarseness of grid-by-number 


(B,) distributions may stem from the B axis having a 


+) 
greater value than their D.- 
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CHAPTER TWO 
THE NUMERICAL METHOD AND THREE EMPIRICAL TESTS 


2.1 General Introduction 

This chapter presents Kellerhals et al (1975) numerical method 
and three empirical experiments which test its predictions for nonuniform 
materials. 

The numerical method, which was explicitly developed for the 
grain-size sampling of thin sections, incorporates solutions to geometric 
and grain-size measure equivalence problems discussed in Chapter One. 
Grid-by-number is taken to be geometrically equivalent to bulk sieve 
analysis (see 1.4.5). Their equations (3a), (3b), and (3c) embody the 
belief that B and D. are systematically related in such a way that 
equivalence is a special case (see 1.4.6 and lore) 

Kellerhals et al (1975) argue that other attempts at determining 
grain size from thin sections of densely packed granular materials have 
run into difficulties for the following reasons: 

1) The sampling procedures employed have not been geometrically 

equivalent to bulk sieve analysis. 

2) Many theoretical approaches employ spherical shaped grains 

or ellipsoids of rotation. Instead, ellipsoids whose k 
values are within the common range should be observed. 

3) Theoretical solutions assume that the grain centers are 

distributed in space according to the Poisson process. 
This is reasonable for dilutely distributed phases but 
not for densely packed granular materials found in 


sediments. The procedures involved in sampling this 
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material have resisted analytical study because there is 
no mathematical definition of a random distribution of 
grain centers, given a closely packed granular phase of 
nonoverlapping grains. 
2.2 The Numerical Method 
2.2.0 Introduction 
The comprehensiveness of the numerical method contrasts with 
the simplicity of the experiments in 1.3. Other than this, any difference 
depends solely upon which is fixed in position, the plane or ellipsoid. 
In 1.3 ellipsoids were moved relative to a fixed plane, whereas, in the 
thin section experiment the fixed plane was transformed into a cutting 
plane operating on a fixed ellipsoid. 
2.2.1 Computation 
In the numerical method, distributions of apparent axes a and 
b were computed by cutting ellipsoids with a large number of planes and 
determining the lengths of the major and minor axes of the elliptical 
intercepts between the planes and the ellipsoid. Cutting planes were 
defined by latitude a, longitude g, and the length r, of the normal from 
the center of the ellipsoid (origin of the spherical coordinate system) 
to the cutting plane (Figure 11). All ellipsoids were assumed to have a 
Square mesh diameter D.=1.0. With Zingg diagram coordinates ke kos 


the relationships between D. and the three axes were given as: 


2 2,. —5 
“f, 201 + 19)] “9, (2) 
Be? (20 ake lard (3) 


C=2k, [2(1 + k5)] 0, (4) 
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Figure 11 


Outline of ellipsoid 
—---- Intersection between ellipsoid 
and coordinate planes 
O Apparent ellipse (intersection between ellipsoid 
and plane) 


<—--> Axes of ellipsoid 


Axonometric Pictorial of an Ellipsoid 
Intersected by a Plane (The true axes 
are A=2C and B=1 .5C A line of sight 
with o=60 and 8=60° is used. The 
plane of intersection is located a 
distance of 0.5C from the center of 
the ellipsoid along a line with o 

=45~ and 8=30°.) (from Kellerhals 
et al 1975, p.83) 
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The procedure produced approximately 8000 to 12000 sample 
ellipses per ellipsoid. The resulting a and b values were arranged in 
20 size classes and plotted as histograms, and as cumulative frequency 
distributions. A total of 49 ellipsoids were sampled using ky and ky 
values 0.15, 0.3, 0.5, 0.67, 0.75, 0.85, and 1.0 (see Kellerhals et al, 
1975, p.83 for further details on numerical method). 


2.2.2 Relationship among De» B and Ky 


Interpretations of equations (2), (3), and (4) requires the 

consideration of the following three points: 

1) D. refers to the side length of a square hole, and not 
the maximum hole size. The maximum dimension of the 
square hole is along its two diagonals whose length 

Lyrv 202 (McGinn, 1971, p.30). 
2) For an ellipsoid the critical dimensions which determine 
whether it can possibly pass through a square hole of a 


given side length Des are its B and C axes. These axes 
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major axis of the 


form an elliptical plane with 9 
parent ellipsoid. 

3) If the elliptical plane formed in 2) is coplanar with the plane 
formed by the sides of the square hole, and the center of 
the square is also that of the ellipse, and the B axis 
of the ellipse falls along one of the diagonals, then three 
cases arise: 
a) The boundary of the ellipse in part or completely 

falls outside the square, thus the parent ellipsoid 


is retained by the square hole (Figure 12a). 


b) The boundary of the ellipse touches the square at 
four points but never falls outside, thus the parent 
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Figure 12 


The Relationship of an 
Ellipse to a Square Hole: 
_ Three Cases 


a) The boundary of the 
ellipse in part or 
completely falls outside 
the square 


b) The boundary of the 
ellipse touches the 
square at four points 
but never falls outside 


c) The boundary of the 
ellipse falls completely 
within that of the 
square 
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ellipsoid is at the threshold of being either retained 
or passed by the square hole (Figure 12b). 


c) The boundary of the ellipse falls completely within 
that of the square, thus the parent ellipsoid passes 
through the square hole (Figure 12c). 

The relationships among D.» ks Weave eandee ine(z), (3o)ecand 
(4) are determined for case b). For illustration, B is calculated using 
(3); where D. is held constant at 1.0 and Ko varied from 0.001 to 1.0. 
As ko increases, the B axis of the ellipse decreases (Table 5 and Figure 
13)..0 At k5=0.001, B=1.4142 which equals La for D.=1.0. For common ko 
values, 0.55-0.75, B:D. ranges from 1.24:1 to 1.13:1 respectively. When 
ko reaches 1.0, B=D.=1.0. This is the special case where the ellipse is 
in fact a circle which touches the square at the midpoints of the four 
Sides. 

2.2.3 Apparent Axial Distributions a and b 

The histograms and cumulative frequency distributions generated 
by the numerical method reveal relationships among ellipsoidal k values, 
true axial values and apparent axial distributions a and b. Spheres and 
ellipsoids of rotation, which have received much attention in the liter- 
ature, have very unusual histograms relative to ellipsoids (k, less than 


1, k, less than 1) (Kellerhals et al, 1975). By means of equations (2), 


Z 
(3), and (4) and Figures 14 a and b which summarize the relationships 
between true axial values, k values and the means cf the apparent axial 
distributions, inequalities can be constructed. For ellipsoids with k 
values within the common range, 0.55-0.75, five pairs of k values were 
selected and the inequalities formed for A, B, C, D,, a, b, and ab 


(mean of a and b). In each case D. was assigned a value of 1.0. The 


inequalities and their numerical values (Table 6) lead to the following 


observations. 
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Table 5 
B Values Calculated by Equation (3), 


D.=1.0, ko range, 0.001 to 1.0 
0. k» B 
1.0 001 1.4142 1.0 
1.0 01 1.4141 
1.0 "02 14139 mee 
1.0 "03 1.4136 
1.0 "04 1.4131 
1.0 05 114124 ah 
1.0 "06 eT 
1.0 07 1.4108 
1.0 ‘08 1.4097 ~ 
1.0 09 1.4085 
1.0 "10 1.4072 
1.0 20 1.3868 aa 
180 "30 1.3546 0. 
1.0 40 1.3131 k, 
1.0 50 1.2649 be 
1.0 60 1.2127 
1.0 "70 1.1586 
1.0 80 1.1043 5 
1.0 "90 1.0512 ; 
a ae 10000 


Te Ose or olemet Ceres) Oe eal ee eS 
B axis length 


Figure 13 
Graph of Table 5 Data 
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Figure 14a Figure 14b 


Difference Between the Sample Mean of Apparent Axes and 
the Corresponding True Axes (from Kellerhals et al 1975, p.88) 
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Table 6 


Axial Inequalities for Ellipsoids 
with Common k Values 0.55-0.75 
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1) Two inequality patterns emerge: 
A>a>B> a ab nC. eb and 
NeooR oweac er fab =aC *) Db 
These inequalities are the same except a and B are reversed 
in position. The reason for this is discernable in Figure 
14a. It reveals that a > B when k, < 0.52 (for all ko) 
and B > a when k, > 0.667 (for all ky). For 0.52 < k, < 
0.667, a > B for lower kK, and higher wo) values, for B>a 
the converse is true. Thus for the five pairs examined 
it is reasonable to expect Ky» i) Hairs: 0.55, 0.59 and 
0.55, 0.75 will exhibit the inequality a > BY 

2) For the five pairs and more generally for any Ki> ko pair 

~with values within the common range, B and a, and C and 

b are closely associated. Within this range 
4.4 > AB2B}100 


> -10.0(Figure 14a) and 


-8.0 -13.8 (Figure 14b). This latter figure 


B-C)100 , 
C 


shows the inequality C > b is true for all pairs of k 
values. 
3) For the five pairs B > Bp 28 C. Inspection of equations 
(2), (3) and (4) demonstrates that B > D. > € whenever 
ellipsoidal k values are less than 1.0 (see also 2.2.2). 
The preceding results were based on grains of a particular 
shape characterized by one point on the Zingg diagram. The sensitivity 
of a distribution to variability in shape was tested by forming a 
combined a-distribution for grains with values clustered around a point 
on the Zingg diagram. This was compared with the distribution of a, for 


grains with the central shape characteristic. Kellerhals et al (1975) 
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found the curves to be very similar and thus concluded that an average 
Shape yalue was justified. 
Zed WInkee Empirical Experiments 

2.3.0 Introduction 

The numerical method simulates the grid-by-number sampling of 
a randomly selected thin section obtained from an isotropic material 
composed of identical ellipsoidal grains. Given the axial dimensions 
of the constituent grains of this uniform material, equations (2), (3), 
and (4) (2.2.1) and Figures 14a and b (25283) predict the D, and mean 
values of apparent axial distributions a and b which would be obtained 
if this material was actually sampled. While results can be extended 
for a material in which grains vary in shape (2.2.3), both grain shape 
and size vary in most clastic sediments. For the numerical method to 
have practical value it must be shown that it is applicable to the 
sampling of these nonuniform materials. 

2.3.1 Experiment 1] 

Kellerhals et al (1975, p.85) compare the median sieve diameter‘ 
predicted by their numerical method with that obtained by sieving, De seyeS0? 
for seven sandstone samples. For all samples, grid-by-number (a,) and 


grid-by-number (b,) distributions were obtained from thin sections. 


) 
Simultaneous grid-by-number and bulk sieve distributions for five of 
these samples were obtained from Friedman (1958, Figure 7). The two 
remaining samples, Y, and Yo were acquired from two artificial sandstone 
blocks made by cementing mixtures of fluvioglacial sand. These blocks 
were then cut to produce thin sections, 

For Friedman's data, bis was converted to C 50° a numerical 


method prediction of C (henceforth: subscript p signifies that the 


term's value depends both on experimental data and numerical method 
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predictions). The conversion of bisg into C 


b we P 
- 150. -. (b-€)100 _ i 
average kp Nes in Figure 14b, C » be C50" (1.0-Y) (bien) (method 


50 is as follows: Using 


used by Kellerhals et al, 1975, p.85) Co the mean counterpart of C50 


can be calculated in the same way except b, and a are used. A similar 


t t 
100b, £5 
but more precise formula for C50 1S TVETOO) ° 


Using C 50 and (4) the predicted median sieve diameter 4.550" 


C p 
p50 2,54 > 
2k, [2(1+k5)] was calculated. The mean counterpart of d.450° d. 


Pp 
can be calculated in the same way except db, ays and Cy are used. The 
qd, 550 of samples Y, and Yo were calculated in a similar manner except 
that (3) was employed. This was made possible because k1> which is 
necessary to the determination of B from Figure 14a), had been estimated 
through microscopic measurement of A. and By axes of 100 constituent 
grains of Y, and Yo. 

The results (Table 7) for the seven samples may be summarized 
as follows. 


1) The k, values of these samples range from 0.59 to 0.69. 


2 
This is within the common range of k values. 


2) at50 > Igp50 > Peso 
ais ” DP aval > D509 except for Cardium Sandstone where 
Ds ieve50 ° 4459 

3) Denso > abeso 
Deievesd > 250 FOr Yo> Osieveso ~ 450 

4) DF eves0 and q5450 values compare reasonably well for the 


seven samples. This may be observed either directly or 


by comparing the ratios 4450°4sp50 and a459:D. 5 aye5q for 


each sample. The mean value of the former is 1.09:1 and 


thet latter racio 1s 1.1021. 
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In general de and D appear to be comparable. Both 


p50 sieve50 
participate analogously in inequalities, form similar ratios and are 


close in actual value. Since q.550 was based on the numerical method 


which assumes a uniform material, and D was produced by sieving 


sieve50 
a nonuniform material, their high degree of correspondence supports the 


hypothesis that the numerical method can be employed effectively in 


Sampling thin section surfaces of nonuniform materials. 


2.3.2 Experiment 2 
Friedman (1958, 1962) discussed the results of thin section 


grid-by-number (a,) and bulk sieve analyses of his 38 sandstone samples. 


Table 8 (adapted from Friedman, 1962, Table 2a) displays a,, D. and 


sieve 


ae 7 ane 
a, 2D. ove for these samples (for Friedman Mean Xs = 709 = fm where 


f = frequency of the different grain-size grades present, and m4 = mid- 


point of each grain-size grade in phi values). 


> * 
For all samples dec VA pee ; 


and actual inequalities in Tables 6 and 7 (sole exception being Cardium 


This conforms with both predicted 


Sandstone where Ap5Q:D =0.96:1). In order to further clarify the 


sieve50 
extent to which the numerical method predicts this inequality to arise, 
at least for uniform materials, and particularly for Kj: kK, pairs within 
the common range, three possible cases (a > De; a=D . and D. > a) were 
investigated. Results (Table 9 and Figure 15a) indicate that the k,> ko 
pairs region where De > a, is very small relative to the a > D. region, 
and is confined to K> Ko pairs of high value. As well, Ky» Ko pairs 
within the common range fall completely within the latter region. 


Assuming that most of Friedman's sandstone samples have mean k values 


within the common range, these results confirm the predictive power of 


*jn 2.3.2 the values of mean terms are in mm. 
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Table 8 
Based on Thin Section Grid-by-Number 


The Ratio a Dee ve 


Analysis and Bulk Sieve Analysis 
(adapted from Friedman, 1962, Table 2a) 


3 Ds ieve ot Dsieve =a 

Sample (@ units) (9 units) (mm) (mm) at sieve 
] 2.55 2.68 .171 E156 reatOet 
2 2310 2.92 154 ~ 132 MA7 21 
5 Fahey) ZeLS LOT. Ay Pa Red 
4 2.39 2.60 191 . 165 eos] 
5 ZN8) 2203 .170 . 141 eet 
6 (dea | 2.99 .147 am WA6) peat We 
7 2205 2401 159 43 deel 
8 20 2.98 .154 Baal lee lel 
9 (esha 7 est 152 . 140 1.09:1 
10 lds eye 2.98 . 142 all ay iPad ay 
1] 0.90 103 550 .490 1.09:1 
12 0.93 1.34 eye) .395 leacoal 
13 0.81 1.24 seyh®) ~423 13521 
14 0.74 eis 599 ~457 tess 
1B) O75 ilinure 595 444 W341 
16 2.50 Lei Sa) OS else 
17 3.64 35/6 .080 .074 10S 1 
18 3.8 3.44 So) .092 Ade 
19 3.02 3.36 mics .097 Ve27e1 
20 Van ote) 2.88 .139 .136 1Te0Ze 
21 Deel LAs AAs .199 1.08: 1 
ZC 2.09 eed 1250 SOLE 1.06:1 
us: (E810 2.4] e203 . 188 1708: | 
24 2.65 2.74 159 .150 120621 
25 2.05 Zei4 » 18241 Aa TO06E1 
26 ea 0 L393 . 308 262 else 
27 3251 SIs .088 [073 Tec] 
28 Saks. Seo ORS . 068 e071 
29 Sey 4.18 .064 7055 1h 1621 
30 3374 3.86 2025 . 069 1.09: 1 
31 3.80 3.94 nOFZ 065 juga ec 
32 Soll Si ee .088 .087 1201 cH 
53 aed 3.42 Pali i: .093 Weeoo 
34 2.74 3703 5 sl22 lege 
35 Zao 2.49 | 57k: 1242) 
36 2205 2.48: 5 ES) .179 1355) 
37 Zac 7209 GAS a | aise 
38 1.88 Zecu eChe .218 12521 


Mean Ratio 1.16:1 
Standard Deviation 0.09 
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Fablen9 te: 
Some Ky ko Pairs for which a:Do=1:1 Or 1.1621 


B from 
‘ky ko Ds equation (3) 
0.667 -1.00 1.00 1.00 
Oeyoy & O82. FAO 1.09 
0.85." 0266. _1.200 1.18 
100 n 0.49) 271200 27 
Ordos l00, +100 1.00 
OP667. "UO 7" #100 He23 
03/57 40243 9471500 12.80 
265700825 12.00 137 


From Figure 14a) 
(a-B)100_ ,  a=B( X_,,) 


B 100 
0 1.00 
="O50 1.00 
-14.9 1.00 
-21.6 1.00 
+16.0 1.16 
- 5.6 1.16 
-10.6 Tal6 
-15.4 1.16 


Figure 15a _ 
Regions where a > D. 
and D. > a 


(adapted from Kellerhals et al 1975, 
Figure 5a) 


Figure 15b 
a:D. Values for kj, kp Pairs 
Within the Conmon Range 
(refer Table 6) 
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the numerical method for nonuniform materials. 

Table 8 indicates that the average 47D. 5 ayer) - 16:1 (standard 
deviation 0.09) for the 38 sample ratios. For the purpose of evaluating 
this ratio with respect to numerical method predictions, the line 
corresponding to a:D.=1.16:1 was drawn (Figures 15a and b). As well, in 
Figure 15b, additional a:D, ratios for common k,> Ko pairs found::in 
Table 6 were plotted. Figure 15b shows that for the a:D. equal width 
"bands" exist which decrease in value (a value decreases) as the kK, and 
Kk» values of a pair approach 1.0. For Ky» ko pairs with common values, 
a:D, ranges from 1.25:1 (k,=k,=0.55) to 1.0331 (k,=k,=0.75). The a:D.= 
1.16:1 line falls centrally within the subregion formed by the common 
range ky» ko pairs. Significantly, the ratios 1.25:1 and 1.07:1, which 


correspond to the EP hey 


either side of the 38 sample mean, also fall within this subregion. 


ratios which are one standard deviation on 


Assuming that the mean k values of most of Friedman's samples fell within 
the common range, these results strongly support the hypothesis that the 
a:D. ratios of the numerical method predict the behavior of ae tave 


for nonuniform materials. 

2.39.0 ‘experiment 3 

In section 1.5 distributions provided by grid-by-number (B,) 
and bulk sieve analysis of exposed gravel bars, were examined. It was 
observed that grid-by-number (B,) distributions were coarser than their 


D counterparts (1.5.1). It was concluded (1.5.4) that there were 


sieve 
two reasons for this inequality arising; geometric non-equivalence due 
to surface paving and the non-equivalence of grain-size measures B and D.- 
These possible causes of non-equivalence can be investigated in detail 


since McGinn's (1971) sampling program differentiates between surface and 
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subsurface deposits and the Kellerhals et al (1975) numerical method 
makes quantitative predictions PneeCninG the relationship among B, De» 
and k» (25252) 4 

McGinn sampled the surface layer with essentially eight different 
procedures, including grid-by-number (B, ) and quadrant (area)-by-weight 
(henceforth: Q is the term denoting the quadrant (area)-by-weight sampling 
procedure whereby all surface grains > 8 mm, within a specified area 
are removed and sieved; Q is the term denoting the mean of the grain- 
size distribution obtained from Q). The subsurface deposit was also 


sampled using bulk sieve analysis (henceforth: D is the term 


Sieve2 
denoting the bulk sieve analysis of the coarse (> 8 mm) portion of the 


subsurface sample; D is the term denoting the mean of the grain- 


sieve2 


size distribution obtained from D ). The weights of Q and D 


sieve2 sieve2 


were then added to provide the D distribution for each sample. 


Sieve 
The advantage of comparing grid-by-number (B,) and Q is that 
only surface grains are being sampled, thus population differences due 
to surface paving are avoided. However, as discussed by Kellerhals and 
Bray (1971) and McGinn (1971), Q is not geometrically equivalent to 
either grid-by-number or volume-by-weight and requires the weighting 
factor 4, (refer to Kellerhals and Bray, 1971, p.1173, 1175 for more 


detail on this procedure; is the term denoting the mean of the 


Q 


sieve 


; grain-size distribution). 
sieve 


and Dejeve? 


sieve’ Dsieve 
both paved and non-paved deposits. If the deposit is paved, it is expected 


The values Q, * can be compared for 


(ahs iL) ae = Da: 
sieve sieve sieve 
In both types of deposits Q > 


9° whereas if non-paved Q SO eteaeiever: 


that sieve 


: : : ] : 
Sieve punSs the weighting factor D applied 


to Q always produces a finer Oe ave distribution, 


*in 2.3.3 the values of mean terms are in mm. 
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Although Q and grid-by-number (B,) are geometrically 


sieye 


equivalent, the numerical method suggests that By > ee More 


specifically, it predicts that for a uniform material composed of 
ellipsoids the B:D. range is 1.4142:1 to 1:1, where B:D.=1.4142:1 as 
ko approaches 0 and B:D.=1:1 for ko=1.0. The B:D. range for common i) 


values. is’ 1.24:1 to 1.13:1, where B:D.=1.24:1 for k,=0.55 and B:D.=1.13:1 


Z 
for k5=0.75(2.2.2). Since the numerical method has been applied quite 


successfully to non-uniform materials in 2.3.1 and 2.3.2 and the grain- 


size measures of Q and grid-by-number (B,) are De and B respectively, 


sieve 
it is expected that the B,: 


: r f i 
etaye range for surface layer gravels will 


be from 1.24:1 to 1.13:1, assuming a common range of ko values, 
The results of McGinn's (1971, Tables 5 and 8) Wilcoxon Tests 


for his 30 samples tend to conform to expectation. It was found that 


(WB) > 0, 


aaa and Q distributions were also 


(grid-by-number (B,) 
similar for median, skewness and kurtosis parameters) and 


Oe =D. 
sieve 


sieve 
(these distributions were also similar for median, skewness and kurtosis 


parameters). It may be inferred that (Q=B, ) > (0.3 jeaD 


sieves The 


inference that By > Ove lends support to the hypothesis that B is a 


coarser grain-size measure than D.. The relative coarseness of Q and 


and D distributions has been predicted above, 


B, compared to Q 


however, their apparent equivalence is not a necessary outcome of any of 


sieve sieve 


these predictions. The equivalence between Q and D distributions 


Sieve Sieve 


implies that the samples may not have been paved. 


Geometric and grain-size measure non-equivalence were further 


examined by calculating the By D and 


; 3 values for McGinn's 
sieve sieve me - 


samples (raw data from McGinn, 1971, Appendix I) and then combining them 


(fablent0) As 


Sieve 


to form ratios By OS i ever? Qsieve sieve2 and Bs! 
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Table 10 cal 
Setiip 2) Merge QD i ever? Be Dsiever 


{sieve sieve? slat BACs eve 
(adapted from McGinn, 1971, Appendix I) 


Sample 
Weight S — ns al $i 

Q: By Qs ieve DSiever a Qsieve’ wabes 
sample Deieye2 (mm) — (mm) (mm ) Deieve2 Ysieve2 Qsteve 
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well, the sample weight ratio Q:D was calculated. If these ratios 


Sieve2 
Show that the surface and subsurface populations were different, even 


though these suspected differences had been compensated for by using 


only grains > 8 mm, then it is clear that Dei aye Was a mixture of two 


distinct populations and the value of D was dependent upon their 


sieve 
grain-size characteristics and proportions. As such, the comparison of a 


D distribution with surface sampling procedure distributions grid- 


sieve 


by-number (B,) or Q » cannot be expected to yield equivalent results, 


sieve 
even if D. and B were equivalent grain-size measures. 


The ratio data in Table 10 clearly indicates that most of 


McGinn's samples were paved and thus D was actually a mixture of 


sieve 
two distinct populations. The surface and subsurface deposits were 


combined in approximately equal proportions since the average value of 


the sample weight ratio Q:D =0.98:1, standard deviation = 0.40. 


sieve2 
The relative coarseness of the surface is confirmed by the average value 


sieve Usieyer | :o3:!> Standard deviation = 0.48. The 


of the ratio 


average values of the ratios B, 2D, =].83:1. standard deviation 0.40, 


sieve2 
and Be Qos aye | 25:1, Standard deviation = 0.25 further emphasize the 
population differences. Grid-by-number (B,) and OT eree which are 


geometrically equivalent surface layer methods, were much closer in 
value than grid-by-number (B,) and De eyer? which were derived from 


surface and subsurface deposits respectively. The Wilcoxon Test 


Reel 
Sieye sieve 


not paved, apparently was a consequence of applying the weighting factor 


result, » which initially suggested that the samples were 


] 1 Ae ‘ é 
D to Q and combining Q with Sieve? values, respectively. At least for 


McGinn's data, these two distinct value-reducing operations performed 


on Q, both produced statistically similar grain-size values. 
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The outcome of the investigation of the ratio B,:Q .. 
t'*sieve for 


McGinn's 30 samples bears heavily on the question of the non-equivalence 
of grain-size measures B and D. since both grid-by-number (B,) and 


Q 


equivalent, but employ different grain-size measures B and Ds. The 


sieve re Surface layer sampling procedures and are geometrically 


sieve 
ko value B:D. range, 1.24:1 to 1.13:1 established by the numerical method, 


average value of the ratio B,: =1.25:1, lies just outside the common 


and corresponds to k=0.53. It is notable that 11 out of the 30 ratio 


values fall within this relatively narrow range, 10 of these being 1.23:1. 


sieve 
deviation 0.25, 1.50:1 to 1:1, corresponds quite closely with the pre- 


The range of the By: values determined by using the standard 


dicted B:D. range, 1.4142 to 1.1. These ratio results uphold the infer- 


ence based on the Wilcoxon Tests that B, >a) implying that B is a 


sieve 
coarser grain-size measure than D.. Furthermore, assuming that the 


mean k values of most of McGinn's samples fell within the common range, 


eave results compare quite well with numerical method predictions. 


a 
The preceding discussion examines why sampling procedure 
experiments which attempted to show the equivalence between bulk-sieve 
analysis and grid-by-number utilizing exposed gravel bars, generally 
found grid-by-number (B, ) distributions to be coarser than those of 


D This effect was probably caused by geometric non-equivalence 


sieve’ 
associated with surface paving and the non-equivalence of grain-size 


measures B and D.. By comparing only geometrically equivalent surface 


layer sampling methods and employing the numerical method, the degree 


sieve 
come confirms the value of the weighting factors associated with the 


to which By > was predicted with reasonable success. This out- 


theory of geometric equivalence (grid-by-number is geometrically 
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equivalent to Q ) and the predictive value of the numerical method's 


sieve 
equation (3) for nonuniform materials. 

2.3.4 Conclusion 

Each of the three sampling experiments test a different facet 
of the numerical method's predictions for nonuniform materials. The 
success of these predictions strongly supports the conclusion that the 


numerical method has broad applicability to the solution of problems 


associated with nonuniform material sampling procedure equivalence. 
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CHAPTER THREE 
TERRACE GRAVEL EXPERIMENT METHODOLOGY 
3.1 General Introduction 

The apparent EGS of a terrace gravel exposure can be treated 
as a thin section surface (1.3.5). However, this surface is different 
from that of a thin section in that both apparent and true axial values 
of grid selected surface grains can be measured easily. Thus terrace 
gravels provide a unique opportunity to test certain numerical method 
predictions not considered in the previous three experiments. 

This chapter presents methods associated with the fourth 
experiment; the grid-by-number sampling of terrace gravels. It 
describes the study area, the field procedures employed in the collec- 
tion of the terrace gravel grid samples, and the methods used to analyze 
the resulting data. 

3.2 The Study Area 

Whitemud and Weed Creeks are north-flowing tributaries of the 
North Saskatchewan River. These two neighbouring stream systems, which 
are found in the vicinity of Edmonton, Alberta (Figure 16), have been 
the subject of geomorphic investigation by Rains (1969) and Shelford 
(1974). Their terrace maps as well as Rains' numerous Whitemud terrace 
stratigraphies, greatly simplified the search for terrace gravels suit- 
able for testing. 

A discussion of the surficial geology of both basins assists 
in clarifying the character of the terrace deposits. In particular, 
their lithological variability is a direct outcome of a complex regional 
surficial geology. The present summary follows Westgate's (1969) 


comprehensive paper on the Quaternary geology of the Edmonton area. 
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The Study Area, Whitemud and 
Weed Creeks near Edmonton, Alberta 
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Bedrock in the Whitemud and Weed Creek drainage basins is 
predominantly Upper Cretaceous in age, and consists of interbedded 
bentonitic shales and sandstones, with some coal seams and bentonite 
beds. Lying unconformably upon the bedrock, and associated with pre- 
glacial valleys, fluviatile preglacial Saskatchewan Sands and Gravels 
are found as terrace and valleyfill deposits. The lithology of this 
material is primarily quartzose sandstone and chert but also includes 
arkosic sandstone, jasper and local bedrock. This composition provides 
strong evidence of Cordilleran origin. Overlying the Saskatchewan Sands 
and Gravels or sitting directly on the bedrock is a lower grayish brown 
till. This and another Laurentide till (upper till) are commonly 
separated by stratified sediments known as Tofield Sand. Both tills 
bear a sizeable proportion of Canadian Shield igneous and metamorphic 
rocks. The northern sections of both basins are veneered with lacustrine 
deposits from the former proglacial Lake Edmonton. 

Shortly after the draining of Lake Edmonton, the North 
Saskatchewan began to cut its valley. Within this postglacial river 
valley four distinct terrace levels can now be observed indicating that 
the river has shown variable rates of degradation. The terrace maps 
of Rains and Shelford revealed analogous terracing patterns in both 
their tributary valleys. The presence of what are termed the "lower, 
middle, upper and higher" cyclic terraces convinced them that periods 
of aggradation and degradation within these tributaries depended upon 
base level changes of the North Saskatchewan, 

The terrace gravels found within the Whitemud and Weed Creek 
basins are lithologically highly variable. The composition is related 


to the local abundance and availability of the Horshoe Canyon Formation, 
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Saskatchewan Sands and Gravels and till-derived materials. Typically, 
the terrace stratigraphy consists of bedrock overlain by a thin stratum 
of alluvial, occasionally imbricated gravel grading upwards into fine 
grained sediments. 

3.3 Field Methods 

3.3.0 Introduction 

This section discusses the factors involved in site selection 

and methods of grid placement, grain selection and measurement. 

S.ae1 “Site Selection 

The terrace exposures examined were those of the lower and 

middle terraces of both Whitemud and Weed Creeks. The following criteria 
were used to determine the site selection. 

1) Site Accessibility - Some promising exposures were 
inaccessible because they were found on quite sheer 
cut-banks. 

2) Dimensions of Gravel Stratum - The exposed terrace 
gravel stratum had to be large enough to accept the 
0.5m square grid, so the exposure's minimum dimension 
had to be at least 0.5M. 

3) Material - Terrace gravel deposits which appeared 
quite isotropic and homogeneous were preferred. In 
most gravel strata anistropy due to imbrication was 
low. Deposits which consisted of substantial 
quantities of shales and sandstones were avoided 
because they were in sheets and were usually highly 


fractured and difficult to extract. 
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4) Surface - The gravel surfaces selected were near 


verti 


initi 


cal and approximately flat. If it was not flat 


ally it was modified so as to decrease photo- 


graphic scale distortion due to distance variations. 


3.3.2 Grid Placement, Grain Selection and Measurement 


In order 
cedure was adopted: 


Step 1) 


Step 2) 


to standardize site investigation, a five step pro- 


A 0.5m grid with 100 intersection points was 
staked firmly in place at the chosen site. Each 
grid was labelled for purposes of photographic 
identification (Figure 17a) e.g., WMM3G; Whitemud 
Creek, middle terrace, site 3, grid G. 

The grid was observed from a fixed frame of 
reference located approximately 1.0m directly in 
front of the grid center (Figure 17b). From this 
vantage point a number of grains were selected 
using the grid intersection points, The approx- 
imate size, shape and location of the selected 
surface grains were recorded on a grid map. 

These systematically gathered grains served as 


the grain population for this grid. 


Two slightly different methods were used to 
select grains. The first method tends to pro- 
duce lower n values (number of grains selected 
per grid) than the second method. It is hypothe- 


sized that these variations in n make no difference. 
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Figure 17b A Grid and Its Fixed Frame of Reference 
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A grain whose minor apparent axis was 

> 8mm (see 1.5.2) and lay under a grid 
intersection point was selected. 

This method required that the grid 
intersection points be examined in a 
specific order. The grid was divided 
into 10 columns each consisting of 10 
grid intersection points. The columns 
were sampled from left to right. Each 
column was examined from the top to the 
bottom. In all cases, if a grain had a 
minor apparent axis > 8mm and lay beneath 
a grid intersection, it was chosen. If 
the grain's minor apparent axis was < 
8mm, the point was considered barren 
and another clast was sought by moving 
either up, right, down or to the left 
along one of the two grid lines composing 
the intersection. The search continued 
as far as the next intersection. Only 
one grain search was allowed per barren 
point, the direction of search being 


varied by 90° clockwise for each successive 


_ barren point encountered. If a suitable 


grain was discovered along the grid line 
and it was not under the adjacent inter- 


section, it was collected as a representative 
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of that intersection. In this latter 
case, as well as when no suitable grain 
was found at all, the barren point was 


left without representation. 


Step 3) Where useful for clarity purposes the surface was 


Step 4) 


Step 5) 


sprayed white (spray paint) and the grid was photo- 
graphed from the frame of reference (Figures 18a and 
b). 

By constant reference to the map the selected grains 
were carefully removed to prevent face collapse and 
subsequent loss. In a number of cases where there 
was a partial face collapse a selected stone could 
be identified and recovered from the debris because 
its size and shape (assistance from the grid map) 
were known and it was partially white due to spraying. 
The A,» B. and C. axes of these selected grains 
were then measured by calipers and recorded. 

Prints of the grid were later examined. Those 
Glass extracted in the field were located by means 
of the grid map and the major and minor apparent 
axes a, and b, were measured. The white paint 


often made the trace outlines of the surface grains 


more evident. 


Since on the prints (approximately 11.5 x 8cm), the 
millimeter divisions could be discerned on the grid 


frame metric rulers, it can be inferred that the 
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Figure 18a Grid and Unsprayed Terrace Gravel Surface 


Figure 18b Grid and Sprayed Terrace Gravel Surface 
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photographic resolution was 1.0mm. Although the 
grids were generally 1/7 life size on the prints, 
and therefore the millimeter divisions approximate- 
Ty 0.14mm on these prints, the calipers could 
measure down to 0.02mm and thus could accurately 
measure the divisions on the millimeter rule. 

3.4 Data Analysis 

3.4.0 Introduction 

The pecrne of this experiment is to observe the effect 
variations in assigning linear dimensions to selected grains have on 
their respective grain-size distributions, and to compare these results 
with those predicted by the numerical method. Because identical grains 
of the surface layer population are employed in a given grid-by-number 
experiment, differences in a grain-size distribution must be a conse- 
quence of treatment (grain-size measurement). 

The remainder of this section discusses the test procedures 
whereby specific predictions of the numerical method are compared with 
the empirical data. Also, a method for examining the effect of sample 
number on the results is described. 

3.4.1 Primary Data Analysis 

The true and apparent axes (A,> By» Cys ay, b,) of each grid- 
selected grain were measured (3.3.2). Four additional grain-size 
measures (AB, » ABC,, BC,, ab, ) (1.5.1) were also computed for each 
grain. The frequency of the data associated with each of these nine 
grain-size measures was determined by number using 0.25 @ class intervals. 
From this grid-by-number data the values for nine cumulative size~- 


frequency distributions were calculated for each grid. 
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The distributions associated with each grid were plotted on 
arithmetic probability paper. For each of these curves, % values were 
determined which permitted the computation of the median, mean, standard 
deviation, skewness and kurtosis parametric values by Folk and Ward 
statistics (King, 1966). 

3.4.2 The Wilcoxon Matched Pairs Signed Ranks Test 

The Wilcoxon Test utilizes both direction and magnitude 
differences within matched pairs to test whether two treatments are 
different. Thigmpained non-parametric test was used because it does 
not require the assumption that the populations being examined are nor- 
mally distributed. 

In the Wilcoxon Test, differences are initially ranked without 
regard to sign. Subsequently, both positive and negative values are 
totalled and the smallest value noted. This value is compared with the 
Wilcoxon Table and its significance determined for the selected probabil- 
ity level (in this experiment a = 0.05 for two-tailed tests; a = 0.025 
for one-tailed tests). The rationale of the test is that if two treat- 
ments are equivalent, then the sums of the positive and negative ranks 
should be about equal. If the sums are considerably different then the 
null hypothesis, Ho» that the treatments do not differ is rejected and 


H,, the alternative hypothesis accepted (see Siegel, 1956, p.75-83 for 


1? 
more details). 

In this study 20 matched pairs were used in each test, each 
pair being derived from a different grid experiment. For a given test 
of two treatment types, all values tested were for the same parameter, 


e.g., the @ values of 20 Beso 4t50 matched pairs. The method of analysis 


is similar to that of McGinn (1971, p.17-20); see also 1.5 and 2.3.3. 
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3.4.3 Testing Procedures A to F 
Test Procedure A 
The Wilcoxon Test (two-tailed) was used to test all possible 
combinations of median and mean matched pairs consisting of both true 
axial and apparent values. The ABC. 59 and ABC, values were matched with 
each of the other five true axial medians and means respectively, and 
tested. As well, the standard deviation, skewness and kurtosis parametric 
values of five pairs (B,-ABC, .B,-a, »ABC,-a, »BC,-ab, .C,-b,) were tested. 
This initial testing procedure serves to detect general 
relationships between true and apparent axial distributions. Special 
emphasis was given to the five pairs because of McGinn's (1971) finding 
that grid-by-number (B,) and (ABC, ) distributions were statistically 
equivalent (1.5.1) and Kellerhals et al (1975, p.84) statement that: 
"Eyen in the case of nonuniform materials, the observed 
distribution of a can be expected to resemble the true 
distribution of B in first approximation, and similarly 
the observed b-distribution will resemble the true C- 
distribution." 
Test Procedure B 
This testing procedure utilized the Wilcoxon Test to examine 
specific numerical method predictions for grid-by-number (B,) and (a,) 
distributions, and grid-by-number (C,) and (b,) distributions. Bisg and 
By were converted by means of Figure 14a) to obtain 4050 and ay values 
respectively; similarly C50 and C, were converted by means of Figure 14b) 


p 
and b. are terms denoting predicted median and mean apparent axial values 


to obtain b 50 and Dy values respectively (henceforth: 450° as bY 50 
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Example of conversion of B50 into a 


+50 : 
average k, = 7—— , average eee 
'@autoo 


p50 


t50 — 
using this average k, and ky in Figure 14a), 22k UE OX 
4550 +50 (ra + by ay is formed in the same way except Ay By and C, 


are used. 


Example of conversion of Ci50 into D560 


C 
average Ky = 7 
t50 (b-C)100 
using this average k, in Figure 14b), ea Ce al Y 


b 550 = Cesolqog * Le Es is formed in the same way except By and C. are 
used. 

The four sets of 20 matched pairs (245072 950° aya, b+ 50750, 
B,-bp) were then compared using the Wilcoson Test (two-tailed). 

Test Procedure C 

This testing procedure utilized the Wilcoxon Test to examine 
the relationship between grid-by-number (C,) and (b,) distributions. In 
contrast to testing procedure B which had similar objectives, bi 50 and b, : 
were converted by means of Figure 14b) to obtain Co50 and a (see 2.3.1] 
for more details on conversion). The two sets of 20 matched pairs 


p50° C.-C ) were then compared using the Wilcoxon Test (two-tailed). 


Test Procedure D 


(C C 


t50 


This testing procedure employed the Wilcoxon Test to investigate 


the relationship among the four methods of estimating the average i) 
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As well, the arithmetic and weighted mean were calculated for 
each of these average ko methods. 

Test Procedure E 

This testing procedure used the Wilcoxon Test (two-tailed) to 
examine the realtionship between predicted sieve median and mean values 
based on either true or apparent axial values. Thus two sets of 20 


matched pairs (0.550 - q. 550° Dep ~ dsp) were compared (henceforth: 


D550 and Us are terms denoting predicted median and mean sieve values 
calculated as follows: 
°+50 


C 
average k, = B we 

t50 B. 

2 


C50 2 ! ~ Ce 2\4 46 : 
p50" Zk, [2(1+k5)] equation (4); DS 2 sts [2(1+k5)] equation (4) 
2 
Calculations required to obtain predicted sieve values d and d 
sp50 b _.SP 
were described in 2.3.1. These required average Ko estimates a 
b t50 
and at respectively. 
ot 


Test Procedure F 

The purpose of this testing procedure is to determine if the 
degree of difference between matched pairs which the numerical method 
predicts are of equal value, is dependent upon either grid sample size 
n, or the use of median versus mean values. 

These particular matched pairs were employed because their 
predicted relationship is one of equivalence (no difference). In 
contrast, the numerical method predicts that differences between matched 
pairs By - a, or BC. 59 ~ bi59 for example, are dependent upon the average 
k values of each grid experiment. In this latter case, differences 
caused by factors such as sample size may be obscured. 


The testing procedure utilized eight sets of 20 matched pairs 
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G C D d 


p50" C50 Rp iep507%sp50.~ “sp50° “+t ~ [p> Mt 
C, - C. and D.. - d_.. The matched pairs were grouped according to the 


size of the sample and by parameter (median or mean) (Table II). 


Table II Matched Pair Groups 


n Median Mean 
20 to 29 group | median group 1 mean 
30 to 39 group 2 median | group 2 mean 
40 to 49 group 3 median group 3 mean 
50 to 59 group 4 median group 4 mean 
60 to 69 group 5 median group 5 mean 


Groups were compared with respect to their minimum and maximum 
differences, range of differences and the weighted mean of their 


differences. 
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CHAPTER FOUR 
TERRACE GRAVEL EXPERIMENT RESULTS AND DISCUSSION 


4.1 General Introduction 

Chapter Four presents the results of the terrace gravel grid- 
by-number experiment described in Chapter Three. The degree to which the 
testing procedures confirm numerical method predictions for this non- 
uniform material is discussed. 
4.2 The Data Used in the Wilcoxon Tests 

The 05, $16, 925, @50, 075, @84,-and G95 values of all true and 
apparent axial distributions are located in Appendix I. With these 9 
values the median, mean, standard deviations, skewness and kurtosis para- 
meters of each distribution can be calculated.* The median and mean 9 
values of all true axial distributions are presented in Appendix II. The 
median and mean @ values of ay and b, apparent axial distributions and 
numerical method predictions can be found in Table 16, Chapter Four. 
4.3 Experimental Results and Discussion 

4.3.0 Introduction 

The results of testing procedures A to F described in 3.4.3 
are presented and examined. 

4.3.1 Testing Procedure A 

The Wilcoxon Test results (Table 12) indicate that only three 
of the true and apparent axial matched pairs tested (B,-ABC, ABC, -a, » 
BC,-ab,) have significantly similar median or mean values. For the five 


matched pairs which were tested using all five parameters it is observed: 


*(Folk and Ward, 1957) 
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Table 12 


Wilcoxon Test Comparison of True and Apparent 


Median 
ABC 59 2459 
A50 HWY 
pees. len in 
B50 Ho 
ABC 50 Hy 
Bisco yy 
Ci50 HOY 
Axial 
Measures 
B, -ABC, 
By -a, 
ABC, -a, 
BC, -ab, 
C,-b, 


Axial Distribution Parameters 


Mean 
ab +50 Peso ABC, = a, aby 
Hy Hy A, H, Hy Hy 
Hy H, AB, Hy Hy H, 
H Hy B. Ho H, H, 
Hy Hy ABC, Ho H, 
Ho H, BC. Hy Hy Hy 
H H C 
] ] t H, Hy H 
Standard 
Deviation Skewness Kurtosis 
Ho Ho Ho 
Ho Ho Hy 
Ho Ho Ho 
Ho Ho Hy 
Ho Ho Ho 


Null Hypothesis. Ho: 


Reject Null Hypothesis. 


the compared axial distributions do not differ 
with respect to the tested grain-size parameter 
at level of significance a=0.05 


H.: the compared axial distributions do differ 
q with respect to the tested grain-size 
parameter at level of significance 
a=0.05 
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B,-ABC, matched pair 

The null hypothesis Ne accepted for all five parameters, thus the two 
distributions are statistically equivalent. These findings are identical 
to those of McGinn (1971). 

Ba, matched pair 

The null hypothesis is accepted for the standard deviation and kurtosis 
parameters. This indicates that while distributions have similar 
characteristics one of the grain-size measures produces coarser results. 
ABC, -a, matched pair | 

The null hypothesis is accepted for all parameters except that of the 
median. While this indicates that the distributions based on these two 
measures are statistically similar the results are slightly confusing 
since B,-ABC, distributions are equivalent. This suggests By-ay distri- 
butions are only slightly different with respect to median and mean 
values. 


BC,-ab, matched pair 


Cat 

The null hypothesis is accepted for all parameters except that of the 
mean. Thus the distributions based on these two different measures are 
quite similar. 

C,-b, matched pair 

The null hypothesis is accepted for standard deviation, skewness and 
kurtosis parameters. Like the By-a, matched pairs the distributions have 
similar characteristics but one of these measures produces coarser 
results. 

Considering that the By-a, distribution median and mean values 


are only slightly different and that the BC,-ab, distributions are quite 


similar, it is probable that the C,-b, distributions have only slightly 
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different median and mean values. 

4.3.2 testing Procedure. B 

The Wilcoxon Tests indicate that matched bole 44507 2p50° 
t50 p50 and B,D, are statistically equivalent. This strongly 
_ confirms the numerical method's predictive capability. 


aap» b 


As discussed in 3.4.3 predicted median and mean apparent axial 
values are derived from median and mean true axial values by using 
Figures 14a) and b). In all cases, these figures predicted that B>a 
and C > b for the average k, and ky values used (Tables 13a and b). Table 
14 which is based on Tables 13a and b, provides information on the ranges 
of X and Y. 
| Since these matched pairs are statistically equivalent it may 


be inferred that the grid-by-number (B.) and (C,) distributions in Testing 


eat 
Procedure A are coarser than the grid-by-number (a,) and (b,) distribu- 
tions, respectively. It appears in approximate terms that the By and 


C, distributions are 10 percent coarser (Table 14). Wilcoxon Tests (one- 


— — 


: : ie Se 

tailed) confirm that Bicg > arso Bt > a, Ci5g > dings and C, > b,- 
4.3.3 Testing Procedure C 

The Wilcoxon Tests reveal that C+50 is equivalent to C50 


but that the null hypothesis must be rejected for Gs and C 


4.3.4 Testing Procedure D 


The Wilcoxon Tests show that of the four average ko es timate 
methods all except — provide statistically equivalent values. The 


arithmetic and wei §hted means support these findings and demonstrate 


|e 


derived ) values are probably slightly greater than the other three 


lee) 


tostimates (Table 15). Wilcoxon Tests (one-tailed) uphold this conclusion. 


* >» here signifies "is coarser than". 
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Table 14 
Ranges of a X and — Y in Tables 13a and b 
Predicted ) 
Term X ¥ Table 
4550 -4,.9 to -16.0 13a 
a, -5.3 to -12.5 13b 
b550 -6.7 to -14.1 13a 
b, -7.7 to -14.1 13b 
Table 15 


Wilcoxon Tests (two-tailed) and Arithmetic and Weighted 
Mean Results for Four Average Ko Estimate Methods | 


Method of estimating C C b b 
average ky and Wilcoxon mar, seo = 20 mie 
Test results By t50 t50 ay 
Arithmetic mean of 

0 ks valucen 0.63 0.59 0.58 0.58 
Weighted mean of 0.63 0.59 0.57 0.58 
20 k, values : 2 : : 


a 


4.3.5 Testing Procedure E 
Wilcoxon Tests indicate that matched pairs 0. 5507 %sp50 and 


D -d._ are statistically equivalent. As discussed in 3.4.3, 0S 550 and 


sp Sp 
D.. values are dependent upon numerical method predictions associated 
sp 
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and d__ are derived 


with true axial median and mean values, while 45550 Sp 


from apparent axial median and mean values by the numerical method. The 
equivalence of these two sets of matched pairs supports numerical method 
predictive capability for nonuniform materials. 


That DS 550 is equivalent to d 550 1s consistent with testing 


C50 _ , Pts0 
procedures B (Creo is equivalent to C 50) and C a and ee 
P t50 t50 
equivalent estimates of the aes ka). 
G C 
; wee Z t50 
Since D550 = [201 + Ng where Sa hG) 
and d = C50 WAG) & Ko)] where k -129 (4), it is reasonable that 
sp50 2k, 2 2 A450 4 
DS 550 is equivalent to q. 550° 
That DS, and Gee are statistically equivalent is also consis- 
tent with testing procedures B (C, is not equivalent to Cy) and 
C, Dy 
C ? — as an estimate of the average ka)» but in a more complex 
ay 
ae 


Using the arithmetic and weighted mean results (testing procedure C) in 


(4); D,,=1.33 (C,), where the average ko=— =0.63 
By ‘t 
=1.41 (C_), where the average k = =, 58. 
p 25 
iv 
Since D, and dee are statistically equivalent then; 


d,,, 


ic 
C, =1.06 Cy 
Thus, for equivalence to occur between pee and de it is probable that 
C, is approximately 6 percent greater than C A Wilcoxon Test (one-tailed) 


supports this conclusion by demonstrating that C, is significantly greater 


hanecee 
than D 
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4.3.6 Testing Procedure F 


The numerical method predicts that the 8 sets of matched pairs 
tested in 4.3 should be equivalent. The Wilcoxon Tests show that 
difference in magnitude and direction for these matched pairs are not 
significant except for C.-C... Table 16 presents testing procedure F's 
n group data and Table 17, an analysis of each group's minimum and 
maximum range, and weighted mean matched pair differences. It is 
observed that: 

1) The maximum difference between matched pairs declines 
from group 1 to group 5. The minimum difference between 
the matched pairs is 0 for all groups. Only for group 
2 mean is the maximum difference greater than that for 

the median based matched pairs. In general, the range 
of differences between matched pairs declines from group 
] to group) 5. 

2) The weighted means of small (.1, .2, .3) and. large (.4, 
.5, .6) differences respectively are fairly similar. 

In the case of group 5 the median and mean weighted 
mean differences are’ less than the others, however, 

this may not be significant since they are based on 

the results of only one grid. The weighted means of 
all difference values (0 to .6) decrease from group 1 

to group 5. Except for group 1, the median difference 
weighted means are greater than their mean counterparts. 

3) Large matched pair differences are limited to a minority 
of grids (Table 18). All groups except groups 3 mean, 


5 median and mean possess at least one grid experiment 
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Table 18 
Grids with Large Differences (.4,.5,.6) 

in Table 17 

Group mee 
group median 0 
group mean ae 
group median 5 
group mean 25 
group median 
group median 
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group median 35 
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with large differences. 

Observation 3) helps to clarify observations 1) and 2) 
concerning maximum differences, maximum range and weighted means (0 to 
.6). Lower n groups have the largest differences, ranges and weighted 
means (0 to .6) because they have a greater proportion of grids with 
large differences and these have the greatest value (i.e., .6). 

In general these results indicate that while small sample 
Size grid tests can be performed such that differences between matched 
pairs predicted to be equivalent are small, there is a tendency for a 
relatively high proportion of these tests to exhibit very large matched 
pair differences. In contrast, a relatively low proportion of the 
larger sample tests tend to display very large matched pair differences. 
These results may be due to a decrease of chance variations with increased 
sample size. 

Observation 2) furnishes evidence that matched pairs derived 
from mean data possess smaller differences. In this case, chance 
variations associated with mean values may be less than medians because 
the means utilize more distribution information. 

4.3.7 Conclusion 

Testing procedures B and C (summarized Table 19) demonstrate 
that the numerical method can be successfully employed to predict 
apparent axial measure values from true axial measure values. The 
converse is also true. Similarly, median and mean sieve diameter 
numerical method predictions based on actual true and apparent grain- 
size values are significantly similar (testing procedure E , Table 19). 
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Table 19 
Summary of Testing Procedures B, C and E 


Tes ting Based on Statistical 
procedure true axial apparent axial relation 
measures measures 


Wilcoxon Test (two-tailed) 


B 4550 4459 no difference 
ay a, no difference 
050 bis no difference 
p by no difference 
c Cts C50 no difference 
ce C, difference 
E D550 q. 550 no difference 
Dsp dsp no difference 
Wilcoxon Test (one-tailed) 
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larger ko values. This latter observation combined with the Wilcoxon 


Test result C, > C, is shown to be consistent with the finding that 
Oe) and d, are statistically similar (testing procedure E, Table 19). 
Generally these results powerfully confirm the usefulness of 
the numerical method as it is applied here to the nonuniform terrace 
gravel deposits. For a given grid experiment the findings of testing 
procedure F indicate that large sample size and the use of mean values 
tends to yield smaller differences between estimates of equivalent 


grain-size measure values. 
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CHAPTER FIVE 
SUMMARY AND CONCLUSIONS 
5.1 General Introduction 

This concluding chapter discusses the numerical method with 
respect to the four empirical experiments examined in this thesis. 
Specific attention is given to the techniques, limitations and advan- 
tages of sampling gravels using only photographs of their apparent 
surfaces. 

5.2 The Four Empirical Experiments and the Numerical Method 

The four empirical experiments, which test numerical method 
predictions for nonuniform materials, vary with respect to material, 
sampling situation and sampling procedure (Table 20). These variations 
permit different aspects of the numerical method to be tested. As 
observed, the results of each of the four experiments confirm numerical 
method predictions. 

The numerical method like any other method, model or theory 
will only achieve popular acceptance as its basic capabilities and 
limitations are clearly defined and thoroughly tested. The four experi- 
ments provide a clear indication of the numerical method's broad predictive 
powers. Duplication or modification of any or all of these experiments 
would further assist in determining the reliability of the numerical 
method predictions. 

Experiments slightly different from the ones described are also 
necessary. Kellerhals et al (1975) predict that if chord length measure- 
ments, such as the maximum chord length in a predetermined direction are 
used instead of apparent axial measurements "...the resulting distribution 


should give directly a close approximation of the D, distribution.". 
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Also Kellerhals et al (1975, p.89) discuss the fact that their 
numerical method predictions have been generated assuming a material is 
isotropic. Since clastic sediments are often anisotropic, the numerical 
method results presented in their paper may not be useful for samples 
with preferred orientation. They propose that this can be corrected for 
by changing the basic mathematical sampling process fundamental to the 
numerical method predictions So as to take into account "...the strength 
of the preferred orientation and the alignment of the thin section(s).”". 
This can be tested. 

5.3 Photographic Sampling of Gravel Surfaces 

5.3.0 Introduction 

Empirical experiments 3 and 4 show that the numerical method 
can be employed successfully in the photographic grain-size sampling of 
both apparent OGS and EGS. This final section reviews the evidence 
supporting this conclusion and discusses the techniques, limitations 
and advantages of sampling each type of surface photographically. 

5.3.1 Photographic Sampling of an Apparent OGS 

McGinn (1971) showed in his exposed gravel bar sampling experi- 
ments (example of an apparent OGS) that grid-by-number (B,) and (b,) 
distributions were quite similar (1.5.2). In Experiment 3 (2.3.3) using 


his data, the numerical method predicted with reasonable success the 


: . It may be inferred from 
sieve 


degree by which B, was greater than 
these results that equal success may be obtained if By was replaced by 

b, Ciners by > ee ale However, as discussed in 1.5.2, the equivalence 
of grid-by-number (B,) and (b,) does depend on the C, axis of most gravel 


bar grains being perpendicular to the surface plane. 
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The numerical method conversion of bisg or b, 


median or mean sieve diameter would be facilitated if an average ko 


to a predicted 


value could be derived from the actual measurement of the true axes of 
50 to 100 grains within the general area where photographic sampling 


is taking place. Using equation (3), the predicted median and mean 


sieve diaméters would equal oe [201 + KS) 2 and t bt + KS) % 
respectively. It must be remembered that these median and mean sieve 
values only apply tc surface layer gravels. 
5.3.2 Photographic Sampling of an Apparent EGS 
In Experiment 4 (Chapters Three and Four) the terrace gravel 
apparent EGS was treated as a thin section surface (equivalence discussed 
in 1.3). In almost every case the median and mean values of grain-size 
measures which were predicted to be equivalent by the numerical method 
were statistically similar (4.3; Table 19). The predicted median and 
mean sieve diameters, ds 950 and a5 can be calculated in the manner 
described in Experiment 1 (2.3.1). Because average ky estimates based 
on either true or apparent axes tend to be comparable (4.3.4, testing 
procedure D) true axes of grains need not be measured in the field. 
Several problems unique to the photographic sampling of terrace 
gravel apparent EGS may be encountered. 
1) The clasts composing the gravel deposit must be 
approximately ellipsoidal in shape. Fluviatile 
quartzite gravel found within the Edmonton area 
is ideal. 
2) The terrace grayel deposit may have a pronounced 


fabric. As discussed in 5.2 numerical method pre- 


dictions may have to be adjusted to take into 
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account preferred orientation. 

3) The area within the grid frame may include more than 
One bed. This should not prove much of a problem as 
long as it is recognized and grain-size distributions 
determined for each bed. 

4) Since much work is done in a cut-bank setting, it is 
often impossible to position oneself properly to take 
the grid photograph. Usually, the photograph is taken 
from approximately one meter from grid center and the 
line of sight passes through this center and is per- 
pendicular to the surface plane. A solution to this 
problem involves placing the camera on a bracket 
attached to the grid frame (Appendix III). 

For terrace gravels, the field and data analysis procedures 

necessary for the determination of the predicted median and mean sieve 


diameters, d and 4, respectively, are provided in Appendix IV. 


sp50 p 
5.3.3 General Features of Apparent Surface Photographic 
Sampling 


There are several points which must be taken into account 
when sampling gravel apparent OGS and apparent EGS. The number of grains 
collected per sample should be quite large, approximately 50 and the 
minor apparent axis of a selected clast should be no smaller than 8nm. 
There may be problems if the gravel is particularly coarse since the 
frame and grid require a relatively flat surface. In this case, a frame 
bearing rulers but no grid may be used. The grid can be superimposed 
later on the photograph and scale determined using the rulers. If grid- 


by-number results are to be directly comparable to those of bulk sieve 
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analysis, the axial values of a clast should be counted as many times 
as it falls under grid intersection points (see Kellerhals and Bray, 
1971, p.1168 and conversion factor theory 1.4.4). 

Ss o44 Advantages of Photographic Sampling 

The grain-size sampling of river bar and terrace gravel 
deposits by photography has a number of advantages. Time savings may 
enable more comprehensive sampling programs to be considered. Because 
of its speed, ephemeral features associated with river channels or 
beaches may be more easily studied. Field work time-savings could be 
paralleled by the automated analysis of the photographs. Finally, it 
may be the only reasonable sampling technique to employ in certain 


undersea or extraterrestial environments. 
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Important @ Values 
of all True and Apparent Axial Distributions 
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Important @ Values 
of all True and Apparent Axial Distributions 
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Important @ Values 
of all True and Apparent Axial Distributions 
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Important 9 Values 
of all True and Apparent Axial Distributions 
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Grid Camera Bracket 

The bracket consists of a shaft which is attached by means of 
a hinge midway along the bottom of the grid frame. When the grid frame 
is in a near vertical position, the shaft is slung by means of two 
equi-length wires connected to the top corners of the frame such that 
the end of the shaft is located approximately over the center of the 
grid. A small camera is placed in a holder which is bolted to the shaft 
near the end. The holder-camera arrangement is approximately one meter 
from grid center and aligned so that the line of sight passes through 
the center and is perpendicular to the surface plane. By means of a 
cable shutter release or timer, this grid frame-camera system may be 
held over terrace gravel deposits, and consistently centered, in-focus, 


photographs taken (Figures 17b, 19). 


Figure 19 Grid Photograph Taken by means of Grid 
| Camera Bracket and Cable Shutter Release 
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Field and Data Analysis Procedures 


Used to Obtain the Predicted Median and Mean 


Sieve Diameters, 4. 550 and deo 


1) Site Selection 


a) 


b) 


c) 


The clasts associated with the terrace gravel deposit 
should be approximately ellipsoidal in shape. 

Because numerical method predictions are for an 
isotropic material, terrace gravel deposits which 
appear quite isotropic and homogeneous are preferable. 
Grain-size analysis should be performed for only one 


bed. 


2) Terrace Gravel Grid Placement and Photography 


a) 


b) 


c) 


d) 


e) 


The gravel surface should be modified so as to be 
approximately flat in order to decrease photographic 
scale distortion due to distance variations. 

ror ease of photography use the grid camera bracket 
described in Appendix III. 

For clarity purposes the gravel surface may be 
sprayed white or colour film used. 

If the terrace gravels are extremely coarse the 
frame with rulers but no grid may be used. The 

grid pattern may be superimposed on the print. 


Each grid should be labelled in some way. 


3) Grain Selection on the Prints 


a) 


50 to 100 grains should be selected. This may require 


more than one grid photograph. 
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b) Select grains = 8mm. found under the grid intersection 
points. If a clast lies under two intersection points 
it must be counted twice (Kel lerhals and Bray, 1971). 
4) Grain Measurement from the Prints 
a) Measure the apparent major axis ay» and apparent 
minor axis bys of each selected grain with calipers. 
b) Determine the scale factor from the grid rulers and 
apply to the apparent axes values. 
5) Primary Data Analysis 
a) The frequency of the apparent major and minor axis 
data respectively is determined by number using 
0.25 9 or 0.50 @ class intervals. 
b) The cumulative size-frequency distribution is plotted 
on arithmetic probability paper. 
c) The @ 16, # 15 and @ 84 values of the distributions 
in b) are noted. The @ 50 values of the apparent 
major and minor distributions are the values of 
4450 and beso? respectively. The Folk and Ward 
mean Pisppare of the apparent major and 
minor distributions yields the mean values a, and 
b, srespectively. 
6) The Numerical Method Median and Mean Sieve Values 
The following presents the method of computing the predicted 
median sieve diameter, d. 550° 


a) Conyert A450 and bi 59 values from @ units to mm. 
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b) Calculate C 


. 
° 


p50 
Des (B - c)100 
Using ky = in Figure 14b, ~——-+—- =_ Y, 
ais C 
C50 (1.0-Y)(b, 50) 
C. is calculated in the same way except Dy and a, are 
used. 
c) Calculate d5550: 
50 2,7 % 
“sp50 2k, [2(1+k5)] 
b 
where k =e 
t50 


dep is calculated in the same except by» a, and Co 


are used. 
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